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— ^  Plane  diffraction  gratings  with  period  2 T lying  in  a  strip 
0  <  y  <  h  in  the  x,y-plane  are  studied.  A  Rayleigh-Bloch  (R-B)  wave  for 
a  grating  is  the  response("^(x,y,p,q)  to  a  plane  wave(  (2Tr)^exp  (i(px  -  qy)  ) 
incident  from  y  >  h  (p  €  R,  q  >  0) .  Thus  (A  +  (p2+  q2))  »  0  in  the 

domain  G  above  the  grating,  Satisfies  the  Dirichlet  or  Neumann  boundary 
condition  on  and  for  y  >  h 
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^+(x,y,p,q)  *  (2ir)"x  exp  (i(px  -  qy)} 
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•/ 


+  l  ct(P.q)  exp  (i  p,  x+  q?  y} 

(p+A)2<p2+q2 


+  I  c^(p,q)  exp  {i  p^x  -  ( Cp  +  if-  P2  -  q2) 1/2  y} 

(p+i.)2>p2+q2 

where  (Pj^q^)  =  (p  +  £>(p2  +  q2  -  (p  +  £)2)^2  )  and  the  summations ^are^ 
over  all  integers  i  satisfying  the  indicated  inequalities^The  paper 
presents  a  construction  of  R-B  waves  and  a  proof  that 
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Introduction. 

The  phenomena  associated  with  the  scattering  of  acoustic  and 
electromagnetic  waves  by  periodic  surfaces  play  a  role  in  many  areas  of 
applied  physics  and  engineering.  Optical  gratings  date  from  the  nineteenth 
century  and  are  still  used  by  spectroscopists .  More  recently,  gratings 
have  been  used  as  coupling  devices  in  integrated  optics  [5].  Trains  of 
surface  waves  on  the  oceans  are  natural  diffraction  gratings  which 
influence  the  scattering  of  electromagnetic  waves  [31]  and  underwater 
sound  [13].  Similarly,  the  surface  of  a  crystal  acts  as  a  diffraction 
grating  for  the  scattering  of  atomic  beams  [14].  The  literature  on 
diffraction  gratings  and  their  applications  is  very  large.  References  to 
work  done  before  1967  may  be  found  in  the  monograph  by  Stroke  in  the 
Handbuch  der  Physik  [29].  More  recent  references  are  given  in  the  papers 
of  Fortuin  [13],  Millar  [20],  Jordan  and  Lang  [16]  and  De  Santo  [9],  among 
others . 

The  first  theoretical  studies  of  scattering  by  diffraction  gratings 
are  due  to  Rayleigh.  His  "Theory  of  Sound"  Volume  2,  2nd  Edition,  published 
in  1896  [22],  contains  an  analysis  of  the  scattering  of  a  monochromatic 
plane  wave  normally  incident  on  a  grating  with  a  sinusoidal  profile.  In 
a  subsequent  paper  [23]  he  extended  the  analysis  to  oblique  incidence. 
Rayleigh  assumed  in  his  work  that  in  the  half-space  above  the  grating  the 
reflected  wave  is  a  superposition  of  the  specularly  reflected  plane  wave,  a 
finite  number  of  secondary  plane  waves  propagating  in  the  directions  of 
the  higher  order  grating  spectra  of  optics,  and  an  infinite  sequence  of 
evanescent  waves  whose  amplitudes  decrease  exponentially  with  distance 
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from  the  grating.  The  validity  of  Rayleigh's  assumption  for  general 
grating  profiles  was  realized  in  the  early  1930's  [11],  following  Bloch's 
work  [4]  on  the  analogous  problem  of  de  Broglie  waves  in  crystals.  Waves 
of  this  type  will  be  called  Rayleigh-Bloch  waves  (R-B  waves  for  brevity) 
in  this  report. 

The  goal  of  Rayleigh's  work  and  the  literature  based  on  it  was 
to  calculate  the  relative  amplitudes  and  phases  of  the  components  of  the 
R-B  waves  and  several  methods  for  doing  this  have  been  developed.  L.  A. 
Weinstein  [32]  and  J.  A.  De  Santo  [6,  7]  have  given  exact  solutions  to 
the  problem  of  the  scattering  of  monochromatic  plane  waves  by  a  comb 
grating;  i.e.,  an  array  of  periodically  spaced  infinitesimally  thin 
parallel  plates  of  finite  depth  mounted  perpendicularly  on  a  plane.  For 
gratings  with  sinusoidal  profiles,  infinite  systems  of  linear  equaltions 
for  the  complex  reflection  coefficients  have  been  given  by  J.  L.  Uretsky 
[30]  and  J.  A.  De  Santo  [8].  Numerical  solutions  of  these  equations  have 
been  given  by  A.  K.  Jordan  and  R.  H.  Lang  [16]  whose  paper  contains 
references  to  numerical  work  by  other  authors. 

The  work  referenced  above  provides  a  satisfactory  understanding 
of  the  scattering  of  the  steady  beams  used  in  classical  spectroscopy. 
However,  modern  applications  of  gratings  in  areas  such  as  integrated 
optics  and  underwater  sound  require  an  understanding  of  how  transient 
electromagnetic  and  acoustic  fields,  such  as  pulsed  laser  beams  and  sonar 
signals,  are  scattered  by  diffraction  gratings.  The  existing  theory  of 
R-B  waves  is  Inadequate  for  the  analysis  of  these  problems. 

The  purpose  of  this  report  is  to  present  an  eigenfunction  expan¬ 
sion  for  diffraction  gratings  in  which  the  eigenfunctions  are  R-B  waves. 
The  theory  can  be  used  to  analyze  the  scattering  of  transient  fields  by 
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diffraction  gratings.  The  analysis,  which  parallels  the  author* s  work 
on  the  scattering  of  transient  sound  waves  by  bounded  obstacles  [34,  35, 
37]  will  be  given  in  a  separate  report. 
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The  theory  of  R-B  wave  expansions  given  below  is  a  generalization 
of  T.  Ikebe’s  theory  of  distorted  plane  wave  expansions  [15],  first 
developed  for  quantum  mechanical  potential  scattering  and  subsequently 
extended  to  a  variety  of  scattering  problems  [2,  19,  25,  26,  27,  35]. 

The  theory  is  based  on  the  study  of  a  linear  operator  A,  called  here  the 
grating  propagator,  which  is  a  selfadjoint  realization  of  the  negative 
Laplacian  acting  in  the  Hilbert  space  of  square  integrable  acoustic  fields. 
The  principal  result  of  this  report  is  a  representation  of  the  spectral 
family  of  A  by  means  of  R-B  waves.  The  R-B  wave  expansions  follow  as  a 
corollary. 

A  key  step  in  developing  R-B  wave  expansions  is  the  introduction 


p.  The  Hilbert  space  theory  of  such  operators  was  initiated  in  a  recent 

article  by  H.  D.  Alber  [3].  Here  Alberfs  powerful  method  of  analytic 

continuation  of  the  resolvent  of  A  is  used  to  construct  the  R-B  wave 

P 

eigenfunctions . 

The  derivation  of  the  R-B  wave  expansions  given  below  is 
restricted,  for  brevity,  to  the  case  of  two-dimensional  wave  propagation. 
Specifically,  the  waves  are  assumed  to  be  solutions  of  the  wave  equation 
in  a  two-dimensional  grating  domain  and  to  satisfy  the  Dirichlet  or 
Neumann  boundary  condition  on  the  grating  profile.  These  problems 
provide  models  for  the  scattering  of  sound  waves  by  acoustically  soft  or 
rigid  gratings  and  of  TE  or  TM  electromagnetic  waves  by  perfectly 


conducting  gratings.  It  will  be  seen  that  the  methods  employed  are  also 
applicable  to  three-dimensional  (and  n-dimensional)  grating  problems. 


Even  with  the  restriction  to  the  two-dimensional  case,  the 
analytical  work  needed  to  derive  and  fully  establish  eigenfunction 
expansions  for  diffraction  gratings  is  necessarily  intricate  and  lengthy. 
This  is  clear  from  examination  of  the  simpler  case  of  scattering  by 
bounded  obstacles  presented  in  the  author's  monograph  [34].  Therefore, 
to  make  the  work  accessible  to  potential  users,  this  report  presents 
only  the  concepts  and  results  of  the  theory,  together  with  the  principal 
ideas  needed  to  derive  them.  Complete  analytical  details  and  proofs  are 
provided  in  a  companion  report. 

The  remainder  of  this  report  is  organized  as  follows.  §1  contains 
the  definitions  of  the  class  of  grating  domains  and  corresponding  grating 
propagators,  §2  contains  the  definition  of  the  R-B  waves  and  their 
classification  into  surface  waves  and  diffracted  plane  waves.  The  concept 
of  the  reduced  grating  propagator  is  introduced  in  §3  and  the  R-B 
surface  and  diffracted  plane  waves  are  shown  to  be  eigenfunctions  and 
generalized  eigenfunctions,  respectively,  of  A^.  The  section  includes 
the  spectral  analysis  of  the  reduced  propagator  A^  corresponding  to  the 
degenerate  grating  whose  profile  is  a  straight  line.  In  subsequent 
sections  the  R-B  wave  expansions  for  general  gratings  are  developed  as 
perturbations  of  this  special  case. 

In  §4  the  analytic  continuation  of  the  resolvent  of  A^  to  a 
suitable  Riemann  surface  is  constructed  by  the  method  of  Alber.  The 
method  leads  to  a  particularly  strong  form  of  the  limiting  absorption 
principle.  In  §5  the  results  of  §4  are  used  to  construct  the  R-B 


eigenfunctions  for  A^  and  to  derive  corresponding  spectral  representations 
and  eigenfunction  expansions  for  Ap.  In  §6  these  results  are  used  to 
construct  the  corresponding  R-B  wave  spectral  representations  and  eigen¬ 
function  expansions  for  the  grating  propagator  A.  §7  contains  concluding 
remarks  concerning  extensions  of  the  theory  and  unsolved  problems. 
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(1.1),  (1.2)  can  be  written 

(1.6)  C  G  C  R2  for  some  h  >  0 
and  the  translation  invariance  (1.3)  takes  the  form 

(1.7)  G  +  (a,0)  -  G 

where  a  >  0  is  the  primitive  period  of  G. 

The  eigenfunction  expansion  theory  for  R-B  waves  that  satisfy 
the  Dirichlet  boundary  condition  is  developed  below  for  arbitrary  grating 
domains.  For  R-B  waves  that  satisfy  the  Neumann  boundary  condition  the 
following  additional  conditions  are  imposed  on  3G,  the  frontier  of  G. 

(1.8)  G  has  the  local  compactness  property,  and 

(1.9)  there  exists  an  x0  G  R  such  that  the  set 
3 G  n  {(x0,y)  |  y  >  0}  is  finite  and  each 
(xQ,y)  in  the  set  has  a  neighborhood  in  R2 
in  which  3G  is  a  regular  curve  of  class  C3. 

Condition  (1.8)  was  introduced  in  [34]  where  it  was  denoted  by  G  G  LC. 

It  is  a  mild  regularity  property  of  3G.  A  simple  sufficient  condition 
for  G  G  LC  is  the  "finite  tiling  condition"  of  [34,  p.  63].  Grating 
domains  that  satisfy  (1.8)  and  (1.9)  will  be  said  to  have  property  S, 
written  G  G  S.  The  class  includes  all  the  piece-wise  smooth  gratings 
that  arise  in  applications.  Examples  Include  the  domains 
G  ■  (X  |  y  >  h(x)}  where  h(x)  is  bounded,  piece-wise  smooth  and  has 
period  a.  A  special  case  is  De  Santo's  comb  grating  for  which 


§1.  Grating  Domains  and  Grating  Propagators. 

The  plane  diffraction  gratings  that  are  studied  in  this  report 
are  the  boundaries  of  the  class,  of  planar  domains  G  defined  by  the 
following  properties. 

(1.1)  G  is  contained  in  a  half-plane. 

(1.2)  G  contains  a  smaller  half-plane. 

(1.3)  G  is  invariant  under  translation  through  a 

distance  a  >  0. 

Domains  with  these  properties  will  be  called  grating  domains.  The  half¬ 
plane  of  (1.2)  is  necessarily  parallel  to  that  of  (1.1)  and  the  transla¬ 
tion  of  (1.3)  is  necessarily  parallel  to  the  edges  of  these  half-planes. 

The  smallest  a  >  0  for  which  (1.3)  holds  is  called  the  primitive  grating 

period.  It  exists  for  all  gratings  except  the  degenerate  grating  for 
which  G  is  a  half-plane. 

It  will  be  convenient  to  introduce  Cartesian  coordinates 

(1.4)  X  -  (x,y)  6  R2 

in  the  plane  of  G  such  that  the  x-axis  is  parallel  to  the  edges  of  the 
half-planes  of  (1.1),  (1.2)  and  to  identify  G  with  the  corresponding 
domain  (open  connected  set)  Gc  R2.  With  this  convention  if 

(1.5)  R2  -  (X  s  R2  |  y  >  c} 

then,  for  a  suitable  orientation  of  the  coordinate  axes,  conditions 
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(1.1),  (1.2)  can  be  written 

(1.6)  R£  C  G  C  R*  for  some  h  >  0 
and  the  translation  invariance  (1.3)  takes  the  form 

(1.7)  G  +  (a,0)  *  G 

where  a  >  0  is  the  primitive  period  of  G. 

The  eigenfunction  expansion  theory  for  R-B  waves  that  satisfy 
the  Dirichlet  boundary  condition  is  developed  below  for  arbitrary  grating 
domains.  For  R-B  waves  that  satisfy  the  Neumann  boundary  condition  the 
following  additional  conditions  are  imposed  on  3G,  the  frontier  of  G. 

(1.8)  G  has  the  local  compactness  property,  and 

(1.9)  there  exists  an  x0  €  R  such  that  the  set 
3G  n  {(x0,y)  |  y  >  0}  is  finite  and  each 
(xQ,y)  in  the  set  has  a  neighborhood  in  R2 
in  which  3G  is  a  regular  curve  of  class  C3 . 

Condition  (1.8)  was  introduced  in  [34]  where  it  was  denoted  by  G  €  LC. 

It  is  a  mild  regularity  property  of  3G.  A  simple  sufficient  condition 
for  G  €  LC  is  the  "finite  tiling  condition"  of  [34,  p.  63].  Grating 
domains  that  satisfy  (1.8)  and  (1.9)  will  be  said  to  have  property  S, 
written  G  €  S.  The  class  includes  all  the  piece-wise  smooth  gratings 
that  arise  in  applications.  Examples  include  the  domains 
G  *  {X  |  y  >  h(x)}  where  h(x)  is  bounded,  piece-wise  smooth  and  has 
period  a.  A  special  case  is  De  Santo ’s  comb  grating  for  which 
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(1.10) 


h(x) 


=  < 


h  >  0  for  x  ■  0 

0  for  -a/ 2  <  x  <  0  and  0  <  x  <  a/2 
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The  Hilbert  space  theory  of  solutions  of  the  wave  equation  in 
arbitrary  domains  G  C  Rn,  developed  by  the  author  in  [33,  34],  provides 
the  foundation  for  the  analysis  of  scattering  by  diffraction  gratings 
given  below.  The  basic  Hilbert  space  of  the  theory  is  the  Lebesgue 
space  L2 (G)  with  scalar  product 

(1.11)  (u,v)  «  [  uOO  v(X)  dX 

JG 

In  addition,  the  definition  of  the  grating  propagators  makes  use  of  the 
Sobolev  spaces 


(1.12)  L“(G)  =  L2(G)  n  (u  |  D®1d“2u  6  L2(G)  for  at  +  a2  <  m}  , 

where  D:  =*  3/3x,  D2  =  3/3y  and  m  is  a  positive  integer,  and  the  space 


(1.13)  l2 (A,G)  =  L*(G)  n  {u  I  Au  6  L2(G)} 

where  A  =  D2  +  D§  is  the  Laplacian  in  R2.  In  these  definitions  the 
differential  operators  are  to  be  interpreted  in  the  distribution-theoretic 
sense  (cf.  [33,  34]). 

The  grating  propagators  for  a  grating  domain  G  are  selfadjoint 

realizations  in  L2 (G)  of  -A,  acting  on  sets  of  functions  that  satisfy 

the  Neumann  or  Dirichlet  boundary  conditions.  These  operators  will  be 
N  D 

denoted  by  A  (G)  and  A  (G) ,  respectively.  Their  domains  are  subsets  of 
L2(A,G)  that  satisfy  the  boundary  conditions  in  a  form  appropriate  to 
arbitrary  domains  G.  In  particular,  functions  u  €  D(A  (G))  are  required 
to  satisfy  the  generalized  Neumann  condition 


for  all  v  G  L* (G)  .  In  fact,  if  one  defines 

(1.15)  L2(A,G)  =  L2(A,G)  O  {u  |  (1.14)  holds  for  all  vG  L2  (G) }  , 

D(AN(G))  *  L2(A,G)  and  AN(G)u  *  -Au  then  AN(G)  is  a  self adjoint 
non-negative  operator  in  L2 (G) .  This  characterization  was  proved  in 

[34].  It  may  also  be  derived  from  T.  Katofs  theory  of  sesquilinear 

/ 

forms  in  Hilbert  space  [17,  Ch.  6].  It  is  known  that  if  3G  is  a  smooth 
curve  then  D(AN(G))  c  L f(G)  and  Vu  has  a  trace  in  L2 (9G)  which  satisfies 
the  Neumann  boundary  condition  [34]. 

To  define  the  grating  propagator  AD(G)  associated  with  the 
Dirichlet  boundary  condition  let 

(1.16)  Lj (G)  =  closure  of  C*(G)  in  L*(G) 
and  define 

(1.17)  L?(A,G)  =  L®(G)  n  L2 (A,G)  , 

D(A°(G))  =*  L2(A,G)  and  AD(G)u  ®  -Au.  Then  Kato*s  theory  of  sesquilinear 
forms  may  be  used  to  show  that  AD(G)  is  also  a  selfadjoint  non-negative 
operator  in  L2 (G) .  Moreover,  it  is  known  that  if  9G  is  a  smooth  curve 
then  every  u  G  (G)  has  a  trace  u | 9G  G  L2 (9G)  and  every  u  €  Lj(G) 
satisfies  u|3G  *  0  [18]. 

N  D 

The  grating  propagators  A  (G)  and  A  (G)  will  be  shown  to  have 

pure  continuous  spectra.  It  follows  that  the  R-B  wave  eigenfunctions 

must  be  generalized  eigenfunctions  which  are  not  in  L2 (G) .  To  define 

N  D 

them  it  will  be  convenient  to  define  extensions  of  A  (G)  and  A  (G)  which 


act  in  the  space 


(1.18)  LzOC(G)  -  P'(G)  n  {u  I  u  6  L2(K  n  G)  for  all  compact  K  C  R2} 

where  O’(G)  is  the  set  of  all  distributions  on  G.  The  following  subsets 
j^o  c 

of  L2  (G)  are  also  needed: 

(1.19)  L?’ZoC(G)  =  L2°C(G)  n  {u  |  D^D^u  €  LjOC(G)  for  at  +  a2  <  m}  , 

(1.20)  lJ’*OC(A,G)  =  L*,£oc(G)  n  (u  |  Au  €  L2oc(G)}  . 

These  linear  spaces  are  all  Frechet  spaces  (locally  convex  topological 
vector  spaces  which  are  metrizable  and  complete  [10])  under  suitable 

J?  Q  Q 

definitions  of  the  topologies.  Thus  L2  (G)  is  a  Frechet  space  with 
family  of  semi-norms 


(1.21) 


PK(U) 


indexed  by  the  compact  sets  KC  R2 .  Similarly,  L®’^OC(G)  is  a  Frechet 
space  with  family  of  semi-norms 


(1.22) 


PK(u) 


|  l  iD^D^2u(X)|2  dX 

^KhG  a14a2<m 

and  L2>^°C(A,G)  is  a  Frechet  space  with  family  of  semi-norms 


1/2 


((  ^  1/2 
(1.23)  p  (u)  -  { | u(X) | 2  +  | Vu(X) | 2  +  |Au(X)|2}  dX 

UKT>G 


The  following  additional  notation  is  used  below; 


(1.24) 


L^om(G)  -  L2 (G)  n  E'(R2) 
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(1.25) 


L‘*Con(G) 


Li(G) 


n 


Lf”(G) 


where  E'(RZ)  denotes  the  set  of  all  distributions  on  R2  with  compact 
supports . 

The  local  grating  propagator  A^*^oc(G)  for  G  and  the  Neumann 
boundary  condition  is  the  extension  of  A^(G)  in  L^°C(G)  defined  by 


(1.26) 


D(A^*^°C(G))  »  lJ,£oc(A,G) 

=  l!,JIoC(A,G)  n  {u 


(1.14)  holds  for  all  v  €  L2,COm(G)} 


and 

(1.27)  AN,£oC(G)u  =  -Au  for  all  u  €  D(AN> 400 (G))  . 


Similarly,  the  local  grating  propagator  A^*^°C(G)  for  G  and  the  Dirichlet 
boundary  condition  is  the  extension  of  AD(G)  in  L2°C(G)  defined  by 

(1.28)  D(AD,AoC(G))  -  L°’Aoc(A,G)  =  L°,Jloc(G)  n  l5’4°C(A,G) 
where 

(1.29)  jD>^oc(g)  =  Closure  of  C~(G)  in  lJ,Aoc(G) 
and 

(1.30)  AD,£oc(G)u  ■  -Au  for  all  u€  D(A0, Aoc(G) )  . 

The  spectral  analysis  and  eigenfunction  expansions  for  A^(G)  and 
A^(G)  are  nearly  identical.  To  emphasize  this,  and  to  simplify  the 
notation,  the  symbol  A  will  be  used  to  denote  either  A^(G)  or  A^(G)  in 

O.OC 

stating  results  that  are  valid  for  both.  Similarly,  the  symbol  A 

will  denote  A^*^°C(G)  or  A^*^oc(G)  except  where  a  distinction  is  necessary. 
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N  D 

The  spectral  theory  of  A  (G)  and  A  (G)  will  be  developed  by 
perturbation  theory,  beginning  with  the  degenerate  grating  R2 *  The 
grating  propagators  for  this  case  will  be  denoted  by 

(1.31)  A*?  =  AN(R2),  A°  -  A°(R*) 

and 

(1.32)  A*>ZoC  -  AN,2'OC(Ris),  aJ’Zoc  -  AD,2,OC(Rj) 

and  the  condensed  notation  A0  for  A^  or  A^  and  Aj0C  for  A^’^oc  or  A^’^oc 
will  be  used. 

The  spectral  analysis  of  A0  can  be  carried  out  by  separation  of 
variables  and  is  essentially  elementary.  Thus  D2  is  essentially  self- 
adjoint  in  L2(R)  with  complete  family  of  generalized  eigenfunctions 
{ ( 2tt )  exp  (ipx)  )  p  E  R}.  Similarly,  D2  and  the  Neumann  boundary 

condition  define  a  selfadjoint  operator  in  L2(0,<»)  with  complete  family 
{(2/tt)1^2  cos  q  y  |  q  >  0},  while  D2  and  the  Dirichlet  boundary  condition 
define  a  second  selfadjoint  operator  in  L2(0,°°)  with  complete  family 
{(2/tt)1^2  sin  q  y  |  q  >  0}.  It  follows  that  the  products 

(1-33)  iJ'oOC.P.q)  “  ^  eip  X  cos  q  y  ,  (p,q)  e  R* 

(1.34)  t(/j(X,p,q)  »  i  eip*x  sin  q  y  ,  (p,q)  €  R* 

are  in  D(A^,^’0C)  and  D(A^’^oc),  respectively,  and  define  complete 

N  D 

families  of  generalized  eigenfunctions  for  A0  and  AQ .  More  precisely, 

N  D 

if  ipQ  is  used  to  denote  either  if/Q  or  ip0  then  the  classical  Plancherel 
theory  can  be  used  to  derive  an  eigenfunction  expansion  and  spectral 
decomposition  for 
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(1.35)  A0  -  I  V  dn0(y) 

J  0 

which  may  be  formulated  as  follows.  First,  for  all  f  €  L2(Ro)  the  limit 

M  M  _ 

(1.36)  f0(p.q)  -  L2(R*)-lim  <P0(X,P,q)  f(X)  dX 

M-*«>  J  o  J  -M 

exists  and 

rM  (-M 

(1.37)  f(X)  *  L2  (Rj)-lim  ip0(X,P,q)  f0(p,q.)  dpdq 

M-*»  J0  J-M 

and 

(1.38)  lf,L2(R?)  "  l^0lL2(Rj)  * 

Moreover,  the  spectral  family  of  A0  is  given  by 

f  -  ■ 

(1.39)  n0(y)  f(X)  -  (|>0(X,p,q)  f0(p,q)  dpdq  . 

{(p.q) |p2+q2<y.q>o) 

Finally,  if  the  linear  operator  $0  :  L2(Rq)  -*•  L2(Rq)  is  defined  by 

A 

$0f  *  f0  then  $0  is  unitary. 

The  principal  result  of  this  report  is  a  generalization  of  this 

eigenfunction  expansion  and  spectral  analysis  that  is  valid  for  the 

D  N 

operator  A  (G)  in  arbitrary  grating  domains  G  and  for  the  operator  A  (G) 

in  grating  domains  G  €  S.  In  these  generalizations  the  R-B  waves  play 

the  role  of  the  eigenfunctions  iJj0  . 


It  will  be  assumed  in  the  remainder  of  the  report  that  the  unit 
of  length  has  been  chosen  to  make  the  grating  period  a  *  2tt.  This 
normalization,  which  simplifies  many  of  the  equations,  does  not  limit 
the  generality  of  the  theory  because  the  general  case  can  be  recovered  by 
a  simple  change  of  units. 

The  definition  of  the  R-B  waves  can  be  motivated  by  considering 
the  reflection  by  a  grating  of  a  plane  wave 

(2.1)  i|>inc(X,p,q)  =  (2-rr)-1  exp  {i(px  -  qy)},  (p,q)  e  R*  . 

inc 

Note  that  the  effect  of  translating  \p  by  the  grating  period  2ir  is  to 
multiply  it  by  a  factor  of  modulus  1: 

(2.2)  *inC(x  +  27T,y,p,q)  *  exp  {2*nrip}  ^inC(x,y ,p,q)  . 

Since  G  is  invariant  under  this  translation  the  reflected  wave,  if  it  is 

inc 

uniquely  determined  by  4/  ,  must  also  have  property  (2.2).  This  suggests 

the 

Definition.  A  function  ip  €  L^,^0C(A,G)  is  said  to  be  an  R-B  wave 
for  G  if  and  only  if  there  exist  numbers  p  €  R  and  co  >  0  such  that 

(2.3)  ip(x  +  2Tr,y)  »  exp  { 2tt ip }  i|Kx,y)  in  G  , 

(2.4)  Aiji  +  •  0  in  G,  and 

(2.5)  \p(X)  is  bounded  in  G. 


If,  in  addition 
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(2.6)  <|>  e  D(AAoc) 

.then  ^  is  said  to  be  an  R-B  wave  for  A. 

The  parameters  u>  and  p  will  be  called  the  frequency  and  x-momentum 
of  the  R-B  wave,  respectively.  Note  that  p  is  only  determined  modulo  1  by 

(2.6) .  The  x-momentum  that  satisfies 

(2.7)  -1/2  <  p  <  1/2 

will  be  called  the  reduced  x-momentum  of  ip*  Property  (2.3)  is  sometimes 
called  quasi-periodicity  or  p-periodicity .  It  is  equivalent  to  the 
property  that 

(2.8)  ^(x,y)  «  exp  {ipx}  <p(x,y)  for  all  (x,y)  e  G 
where 

(2.9)  <p(x  +  27T,y)  *  <P(x,y)  for  all  (x,y)  €  G  . 

Solutions  of  the  Helmholtz  equation  (2.4)  are  known  to  be 

analytic  functions.  In  particular,  each  R-B  wave  for  A  satisfies 
^  6  C°°(G).  Hence,  the  function  $  in  (2.8)  is  in  C°°(R^)  and  has  period 
27 t  in  x.  In  follows  from  classical  convergence  theory  for  Fourier  series 
that  ip  has  an  expansion 

(2.10)  ip(x, y)  -  l  ipf(y)  exp  {i(p  +  £)x),  (x,y)  €  R*  , 

tez  4  n 

where  2  denotes  the  set  of  all  integers.  The  series  converges  absolutely 
and  uniformly  on  compact  subsets  of  R^.  Moreover,  the  partial  derivatives 
of  \p  have  expansions  of  the  same  form  which  may  be  calculated  from  (2.10) 
by  term-by-term  differentiation  and  which  have  the  same  convergence 


properties.  It  follows  that  the  coefficients  t/^(y)  in  (2.10)  must 
satisfy 
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(2.11)  4»£(y)  +  (u>2  -  (p  +  £)2)  <|^(y)  *  o  for  y  >  h  . 

Hence  the  terms  in  the  expansion  (2.10)  have  the  following  forms, 
depending  on  the  relative  magnitudes  of  to  and  |p  +  &| . 

to  >  1  p  +  &1  .  In  this  case  there  exist  constants  and  c^  such 

that 

(2.12)  ^(y)  exp  {i(p+ £)x}  **  c^exp  {i(p^x+ q^y)  }  +  exp  {i(p^x  -  q^y)  } 
where 

(2.13)  p^  =  p  +  £,  q&  =  (id2  -  (p  +  £)2)1/2  >  0  . 

the  two  terms  in  (2.12)  describe  plane  waves  propagating  in  the  directions 
(P£ >*<?£)<  Since  +  *  to2  these  vectors  lie  on  the  circle  of  radius  to 

with  center  at  the  origin  and  their  x-components  differ  by  integers. 
Clearly  there  are  only  finitely  many  such  terms. 

0)  <  Ip  +  &1  .  In  this  case  i/^(y)  is  a  linear  combination  of  real 
exponentials  in  y  and  the  boundedness  condition  (2.5)  implies  that 

(2.14)  ^(y)  exp  {i(p+  A)x>  -  c^  exp  {-((p  +  l)2  -  to2)  1/2  y}  exp  (i(p  +  )L)x} 

where  ((p  +  £,)2  -  oo2)1^2  >  0.  In  the  application  to  diffraction  gratings 
terms  of  this  type  will  be  interpreted  as  surface  waves. 

to  ■  1  p  +  i\  .  In  this  limiting  case  ^(y)  is  a  linear  combination 
of  1  and  y  and  (2.5)  implies  that 

(2.15)  ^(y)  exp  “  ci  exP  fi(P  +  £)x)  • 


18 


Physically,  (2,15)  describes  a  plane  wave  that  propagates  parallel  to 
the  grating;  i.e.,  a  grazing  wave.  These  waves  divide  the  plane  waves 
(2,12)  from  the  surface  waves  (2,14).  The  frequencies  {u)»  |p+Jl|  |  Z} 
are  called  the  cut-off  frequencies  for  R-B  waves  with  x-momentum  p. 

An  R-B  wave  ij;  for  G  (for  A)  which  satisfies  the  additional 
conditions 

(2.16)  c^  *  0  (resp.  c£  *  0)  for  all  l  such  that  a)  >  |p  +  l\ 

will  be  said  to  be  an  outgoing  (resp.,  incoming)  R-B  wave  for  G  (for  A). 
If 

(2.17)  c£  ■  c+  »  0  for  all  l  such  that  u)  >  |p  +  Jl| 

then  ^  will  be  said  to  be  an  R-B  surface  wave  for  G  (for  A) .  Of  course 

an  R-B  surface  wave  for  A  is  both  an  outgoing  and  an  incoming  R-B  wave 

for  A.  It  is  interesting  that  these  are  the  only  outgoing  or  incoming 
R-B  waves  for  A,  This  is  a  consequence  of 

Theorem  2.1.  Every  outgoing  (resp.,  incoming)  R-B  wave  for  A  is 
an  R-B  surface  wave  for  A. 

A  proof  of  this  result  has  been  given  by  Alber  [3]  in  the  case 
where  3 G  is  a  curve  of  class  C2 .  The  method  is  to  apply  Green’s  theorem 
to  the  R-B  wave  for  A  and  its  conjugate  in  the  region 
G  n  {X  |  -tt  <  x  <  7T,  y  <  R).  In  the  case  of  an  outgoing  R-B  wave  for  A 

this  yields  the  equation 

(2.18)  l  (w2  -  (p  +  A)2)1/2  |c+|2  -  0 

co>  |  p+Jt  |  K 

which  implies  that  *  0  when  uj  >  |p  +  l\ .  For  general  grating  domains 
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the  application  of  Green's  theorem  must  be  based  on  the  generalized 
boundary  conditions,  as  in  [34,  p.  57]. 

It  will  be  seen  in  §4  that  diffraction  gratings  may  indeed 
support  R-B  surface  waves  and  the  question  arises  whether  geometric 
criteria  for  the  non-existence  of  such  waves  can  be  found.  In  the  case 
of  the  Dirichlet  boundary  condition  such  a  criterion  was  found  by  Alber 
[3]  by  adapting  a  method  of  F.  Rellich  [24]  and  D.  M.  Eidus  [12].  Spe¬ 
cialized  to  the  grating  domains  considered  here,  Alber 's  theorem  implies 

Theorem  2.2.  Let 


(2.19) 


G  =  {X  |  y  >  h(x)  for  all  x  €  R} 


where  h  G  C2  (R)  and  h(x  +  2tt)  »  h(x)  for  all  x  G  R.  Then  AU(G)  has  no 
R-B  surface  waves. 

Theorem  2.1  implies  that  R-B  waves  for  A  may  be  determined, 
modulo  R-B  surface  waves,  by  specifying  either  the  coefficients  c^  with 
a)  >  |p  +  £|  (the  incoming  plane  waves)  or  the  coefficients  c"£  with 
a)  >  |p  +  l\  (the  outgoing  plane  waves).  R-B  waves  for  A  that  contain  a 
single  incoming  or  outgoing  plane  wave  will  be  used  in  the  R-B  wave 
expansions  given  in  §6  below.  These  are  the  grating  waves  originally 
introduced  by  Rayleigh.  Physically,  they  are  the  wave  fields  produced 
when  the  grating  is  illuminated  by  a  single  plane  wave.  Here  they  will 
be  called  R-B  diffracted  plane  wave  eigenfunctions  for  A  or,  for  brevity, 
R-B  wave  eigenfunctions  for  A.  There  are  two  families  determined  by  the 
presence  of  a  single  incoming  and  outgoing  plane  wave,  respectively. 

The  plane  waves  4>inc(X,p,q)  and  ipinc(X,p , -q)  defined  by  (2.1)  are 
incoming  and  outgoing  R-B  waves,  respectively,  with  x-momentum  p  and 
frequency 


. . . . . 1 
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(2.20)  .  u  -  U>(p,q)  -  (p2  +  q2)l/2  . 

The  scattering  of  these  waves  by  a  grating  will  produce  outgoing 
(resp.,  incoming)  R-B  waves  with  the  same  x-momentum  and  frequency. 

Hence  the  R-B  wave  eigenfunctions  may  be  defined  as  follows. 

Definition.  The  outgoing  R-B  diffracted  plane  wave  for  A  with 
momentum  (p,q)  €  R*  is  a  function  i|>+(X,p,q)  such  that 

(2.21)  ^+(*,P,q)  is  an  R-B  wave  for  A,  and 

(2.22)  ^+(X,p,q)  -  ^inC(X,p,q)  +  ip®C(X,p,q) 
sc 

where  is  an  outgoing  R-B  wave  for  G.  Similarly,  the  incoming  R-B 
diffracted  plane  wave  for  A  with  momentum  (p,q)  €  R*  is  a  function 
i|;_(X,p,q)  such  that 

(2.23)  'J,_(*  »P»q)  is  an  R-B  wave  for  A,  and 

(2.24)  l|>_(X,p,q)  =  <|>inC(X,p,-q)  +  i|>fc(X,p,q) 

where  is  an  incoming  R-B  wave  for  G. 

The  uniqueness  of  ip+(X,p,q)  modulo  R-B  surface  waves  follows 
from  Theorem  2.1,  as  was  remarked  above.  Their  existence  for  the  class 
of  gratings  defined  in  §1  is  proved  in  §6.  Note  also  that  the  defining 
properties  imply  that 

(2.25)  ^_(X,p,q)  -  ^+(X,-p,q)  . 

Hence  the  existence  of  the  family  \p  follows  from  that  of 

In  the  half-plane  above  the  grating  the  R-B  waves  have 
Fourier  expansions  (2. 10) .  For  the  function  ^  the  expansion  has  the  form 
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^+(x,y,p,q)  =  ( 2tt) ~ 1  exp  {i(px  -  qy)  } 

(2.26)  +  l  ct(p,q)  exp  {i(p?x  +  q*y) } 

(p+Jt)2<p2+q2 

+  [  c^(p,q)  exp  {ip^x}  exp  {-((p+ X-)2  -  p2  -  q2)  1/2  y} 

(p+fc)2>p2+q2 

where 

(2.27)  (pt,q4)  -  (P  +  Mp2  +  q2  “  (P  +  ^)2}l/2  )  €  R2 
defines  the  momentum  of  the  reflected  plane  wave  of  order  Z.  Similarly, 

i|>_(x,y,p,q)  =  (2tt)-1  exp  (i(px  +  qy) } 


(2.28)  +  l  c^(p,q)  exp  {i^x  -  q£y)  >  ; 

(p+Jl)2<p2+q2  j 

+  1  c^(p,q)  exp  { ip^x}  exp  (-((p+ £)2  -  p2  -  q2) 1/2  y) 

(p+Jl)2>p2+q2 

The  relation  (2.25)  implies  that  the  coefficients  c^(p,q)  in  (2.26),  * 

(2.28)  satisfy  | 

I 

-  i 

(2.29)  c~(p,q)  »  c^(-p,q)  for  all  (p,q)  €  R2  and  Z  6  Z  .  \ 

i 

i 

The  surface  wave  terms  in  (2.26)  and  (2.28)  are  exponentially  j 

decreasing  functions  of  y  except  when  the  wave  frequency  ! 

w(p,q)  *  (p2  +  q2)1^2  *  |p  +  H|  for  some  Z  €  Z.  These  are  precisely  ! 

the  cut-off  frequencies  mentioned  above.  In  momentum  space  they  form  • 


» 


the  exceptional  set 


E  is  a  set  of  confocal  parabolas  with  foci  at  (0,0),  axes  along  the 
p-axis  and  directrices  p  +  A  *  0,  A  €  Z.  Two  members  of  the  family  with 
directrices  p  +  A  =  0,  p  +  m  *  0  are  disjoint  if  A  and  m  have  the  same 
sign  and  intersect  orthogonally  if  A  and  m  have  opposite  signs*  The 
family  E  thus  divides  R*  into  a  system  of  curvilinear  rectangles. 

In  the  special  case  of  the  degenerate  grating  Rq  comparison  of 
(1.33),  (1.34)  with  (2.26),  (2.28)  shows  that  for  the  Neumann  case 
=  i|>+  =  c“(p,q)  =  (2tt)  1  and  all  other  c^(p,q)  «  0.  Similarly, 

for  the  Dirichlet  case  =  *  -i^  ,  Co(p,q)  =  -(2tt)~1  and  all  others 

^(P*^)  =  0.  Thus  in  these  cases  there  is  no  scattering  into  higher 
order  grating  modes  or  surface  waves,  as  was  to  be  expected.  Note  that 
the  defining  properties  (2.22),  (2.24)  can  be  rewritten  as 

(2.31)  MX.p.q)  -  ^0(X,p,q)  +  l|^(X,p,q) 

where  is  defined  as  at  the  end  of  §1  and  ^  and  are,  respectively, 
outgoing  and  incoming  R-B  waves  for  G.  This  decomposition  exhibits  the 
R-B  wave  eigenfunctions  for  G  as  perturbations  of  those  for  Rq .  The 
decomposition  is  used  below  for  the  construction  of  and  the  derivation 
of  the  eigenfunction  expansions. 


§3 .  The  Reduced  Grating  Propagator  Ap . 

The  quasi-periodicity  property  (2.3)  of  the  R-B  waves  implies 
that  they  are  completely  determined  by  their  values  in  the  domain 

(3.1)  n  =  G  n  (X  I  -IT  <  X  <  It}  . 

Moreover,  (2.3)  and  the  equation  obtained  from  it  by  x-dif f erentiation 
define  boundary  conditions  that  must  be  satisfied  by  R-B  waves  on  the 
portions  of  dft  where  x  =  ±TT.  These  observations  are  used  below  to  show 
that  the  R-B  surface  waves  and  diffracted  plane  waves  for  G  are  eigen¬ 
functions  and  generalized  eigenfunctions,  respectively,  of  a  p-dependent 
self ad joint  realization  of  -A  in  L2 (ft) .  This  operator,  which  will  be 
denoted  by  A^  and  called  the  reduced  grating  propagator,  provides  a 
basis  for  the  construction  of  the  R-B  waves  for  G. 

The  definition  of  the  grating  domains  in  §1  implies  that  the 
reduced  grating  domains  Q  satisfy 

(3.2)  Bh  C  Q  C  Bq  for  some  h  >  0 
where 

(3.3)  Bc  a  n  {X  |  -7T  <  x  <  7t)  *  {X  |  -7T  <  x  <  ir,  y  >  c}  . 

The  notation 

(3.4)  Y  *  {y  |  ("TT , y)  6  G)  =  (y  I  (-ff.y)  e  G} 

will  also  be  used.  The  definition  of  the  reduced  grating  propagators 

M  T) 

Ap(Q)  and  A^(ft)  associated  with  Q  and  the  two  boundary  conditions  will 
be  based  on  the  function  space 
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(3.5)  I.2,P(n)  =  L*(ft)  n  (u  |  u(ir,y)  »  exp  {2irip}  u(-ir,y),  y  e  y) 


Sobolev's  imbedding  theorem  [1]  implies  that  every  u  €  Lj (fi)  has  boundary 

values  u(±ir,y)  in  Lz°C(y)  and  L*,P((2)  is  a  closed  subspace  of  L*(ft). 

N 

The  operator  A^(ft)  is  defined  by 
D(A^(ft))  =  L*’P(£2)  O  Li  (A,fi) 


(3.6) 


n  {u 


ft 


{ (Au)v+  Vu  •  Vv}  dX  =  0  for  v  €  L*  ,P(ft) } 


N  N 

and  A  (ft) u  *=  -Au.  It  can  be  shown  that  A  (ft)  is  the  selfadjoint  non- 
P  P 

negative  operator  in  L2 (ft)  associated  via  Kato’s  theory  with  the  sesqui- 
linear  form  defined  by  the  Dirichlet  integral  acting  on  the  domain 
L2,^(ft).  By  applying  elliptic  regularity  theory  [1]  and  Sobolev’s 
imbedding  theorem  it  can  also  be  shown  that  every  u  €  D(Ap(ft))  satisfies 
the  p-periodic  boundary  conditions 

u(TT,y)  *  exp  { 2tt ip }  u(--rr,y),  y€  y 


(3.7) 


D1u(iT,y)  *  exp  {2iTip}  D1u(-iT,y),  y  €  y 


Moreover,  if  3G  is  a  smooth  curve  then  it  follows  from  (1.14)  as  in  §1 

N 

that  functions  u  €  D(A^(ft))  satisfy  the  Neumann  boundary  condition  on 


(3.8) 


r  a  3G  Pi  ft 


where  ft  is  the  closure  of  ft  in  R2 . 

To  define  A^(ft)  several  additional  function  spaces  are  needed. 
The  subset  of  C°°(G)  consisting  of  functions  that  satisfy 

(3.9)  <(>(x  +  27r,y)  *  exp  {27Tip}  <J>(x,y)  for  all  (x,y)  €  G 
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(3.10)  supp  <f>  C  G  n  {X  |  y  <  p}  where  p  =  p (<J>) ,  and 

(3.11)  dist  (supp  0 , 3G)  >  0  , 
will  be  denoted  by  c”(G) : 

(3.12)  C*(G)  =■  C°I>(G)  n  {0  |  (3.9),  (3.10)  and  (3.11)  hold}  . 

P 

The  restrictions  of  such  functions  to  0  defines 

(3.13)  C"(fi)  -  -  $\Q  |  4>  €  C”(G)} 

Finally, 

(3.14)  L2,P(J))  =  Closure  in  L*(n)  of  c”(ft)  . 

The  operator  Ap(ft)  is  defined  by 

D(A°(fi)) 

(3.15) 

-  L?,P(fi)  n  L^(A,fi)  n  {u  |  [  { (Au)v+  Vu  •  Vv}  dX  *  0  for  v€  ,p(ft)  1 

Jn 

and  A^(0)u  =  -Au.  In  this  case  it  can  be  shown  that  A^(ft)  is  the  self- 
P  P 

adjoint  non-negative  operator  associated  via  Kato's  theory  with  the 

sesquilinear  form  defined  by  the  Dirichlet  integral  acting  on  the  domain 

Again,  functions  in  D(A^(ft))  satisfy  the  p-periodic  boundary 

conditions  (3.7).  Moreover,  if  8G  is  a  smooth  curve  then  functions 

u  £  D(A^(ft))  satisfy  the  Dirichlet  boundary  condition  on  T. 

N  N 

The  operators  A^(ft)  and  A^(Q)  will  be  shown  in  §4  to  have 
continuous  spectra.  To  define  corresponding  generalized  eigenfunctions 
it  will  be  convenient  to  define  extensions  of  A^(ft)  and  A^(ft)  in  (ft)  • 

£oc 

The  following  subsets  of  L2  (Q)  are  also  needed: 


(3.16)  l2,P>  °C(ft)  *  Lz*  n  (u  |  u(7r,y)*exp  {2nip}  u(-7T,y),  y€  y} 

(3.17)  L2,P,A°C(fi)  -  Closure  of  C*(S1)  in  L*,£oC(n)  . 

j  5»oc 

Each  is  a  closed  subspace  of  the  Frechet  space  L2 9  (Q) .  The  sets 

[  L;’p’com(fi)  =  L^P(n)  n  L~m(n) 


(3.18) 


.  D,p, com 


m 


L?,p(n)  n  L^om((2) 


will  also  be  used. 

The  operator  Ap,^oc(fi)  is  the  extension  of  A^Cft)  in  L2°C(^) 
defined  by 


_  ,  AN,£oc/m.  Tl,p,£oc,ON  _  i,i,oc/A  nN 
D(A  *  (ft))  =  L2  (ft)  n  l2’  (A, Si) 


(3.19) 


n 


(u  |  [  {  (Au)v  +  Vu  •  Vv}  dX  *  0  for  v€  L 1  ,P  ,com< 

'ft 


and  AN,^°C(Si)u  *  -Au.  Similarly,  A^,^'0C(Q)  is  the  extension  of  A^(ft)  in 
P  P  P 


L2°C(ft)  defined  by 


D (A^ ' ^°C (£1) )  -  n  L-l0<:(4, 


fl) 


(3.20) 


n  {u 


{  (Au)v+  Vu  •  Vv}dX«0  for  v  €  L?’P,COffl(n)  : 

a 


and  AD,£oc(Q)u  *  -Au.  It  is  easy  to  verify  that  D(A^’£oc(Q))  and 
P  P 

D(A^>£oc(fl))  are  closed  linear  subspaces  of  the  Frechet  space  L1,£oc(A,$ 
P  2 

and  hence  are  themselves  Frechet  spaces . 

The  reduced  grating  propagators  for  the  degenerate  grating  will 


be  denoted  by 
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(3.21) 

(3.22) 


AN  -  AN(Bft),  A°  -  AD(B0),  and 

0,P  p  0 '  *  0  ,p  p  o 


^,loc  _  aN,£oc 
'  0  .P 


A*’  — (B0),  A0>p 


D,2oc  D,2oc 


(B0 )  • 


Moreover,  the  condensed  notation  of  §1  will  be  used;  i.e.,  Ap  will  be 

used  to  denote  either  Ap(ft)  or  A  (ft)  in  stating  results  valid  for  both. 

Similarly,  Ap0C  will  be  used  to  denote  Ap,^'0C(ft)  or  Ap’^0C(ft).  In 

particular,  for  the  degenerate  grating  the  notation  Aq  is  used  for 

Aa  and  Aa  and  A  is  used  for  A^ ’  and  A^  _ 

0 ,P  o»P  o,p  o»P  0,p 

Note  that  all  the  p-dependent  function  spaces  defined  above  are 
periodic  functions  of  p  with  period  1.  It  follows  that 

,£oc  Aioc 


(3.23) 


A  A  ,  AZIZ  -  A - for  all  m  €  Z  . 

p-hn  p 9  p+m  p 


£oc 

Hence  it  will  suffice  to  study  A  and  A  for  the  reduced  momenta 

7  V  P 

p  €  (-1/2, 1/2] . 

The  resolvent  set  and  spectrum  of  A  will  be  denoted  by  p(A  ) 

p  P 

and  a(Ap)  respectively.  Clearly  a(Ap)  C  [O,00)  since  Ap  is  selfadjoint 
and  non-negative.  In  fact,  it  will  be  shown  that 


(3.24)  a(Ap)  «  [p2,«)  for  all  p  e  (-1/2, 1/2]  . 

This  was  proved  directly  by  Alber  in  the  cases  considered  by  him  [3]. 
Here  it  follows  from  the  eigenfunction  expansions  for  Ap  given  in  §5. 
a(Ap)  is  a  continuous  spectrum  which,  in  general,  will  have  embedded 
eigenvalues.  It  will  be  shown  In  §4  that  a0 (Ap) ,  the  point  spectrum  of 
Ap,  is  discrete;  that  is,  each  interval  contains  finitely  many  eigen¬ 
values  of  Ap  and  the  eigenvalues  have  finite  multiplicity.  It  is  of 
interest  for  the  applications  to  diffraction  gratings  to  have  criteria 
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for  o0( Ap)  to  be  empty.  While  completely  general  criteria  are  not  known 
it  will  be  shown  that  the  hypotheses  of  Theorem  2.2  imply  0o(Ap)  *  4> 
for  all  p  €  (-1/2,1/21. 

Eigenfunction  expansions  for  are  derived  in  §5  by  perturbation 
theory  starting  from  Aq  .  The  expansions  for  Afl  p,  which  are 
elementary,  are  recorded  here  as  a  starting  point  for  the  analysis  of  Ap. 
Separation  of  variables  applied  to  A^  leads  to  the  complete  family  of 
generalized  eigenfunctions 

(3.25)  <j)j  (X,p+m,q)  =  j  ei(P+m)  x  cos  qy,  m  €  Z,  q  >  c 

where  p  e  (-1/2, 1/2]  is  fixed.  Similarly,  for  one  finds  the  complete 

family 

(3.26)  <J>o  (X,p+m,q)  -  ^  ei^p+m^*X  sin  qy,  m  e  Z,  q  >  0  . 


To  describe  the  eigenfunction  expansions  for  A  the  condensed  notation 

o  ,p 

N  D 

(j)0  (X,p+m,q)  will  be  used  to  denote  either  <)>0  or  <pg .  Note  that 

(3.27)  <j>0(X,p+m,q)  *  (X,p+m,q)  |  , 

that  is,  the  generalized  eigenfunctions  for  A  are  obtained  from  those 
of  A0  by  restricting  X  to  B0  and  the  x-momentum  parameter  to  the  lines 
p1  »  p  +  m  with  m  e  Z  and  p  e  (-1/2, 1/2]  fixed.  Classical  Plancherel 
theory  implies  that  if  R0  ■  (0,»)  then  for  all  f  €  L2(B0)  the  limits 


f0(p+m,q) 


L2(R0)-lim 


B 


o  » 


(3.28) 


$0 (X,p+m,q)  f(X)  dX 
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exist  for  m  e  Z  and  p  €  (-1/2, 1/2]  fixed,  where  B0  M  *  B0  n  {X  |  y  <  M). 
Note  that  the  L2 (RQ) -convergence  refers  to  the  variable  q.  Moreover, 
Parseval's  formula  holds  in  the  form 


(3.29) 


l 

mEZ 


|f0(p4m,*>l[2(Ro)  • 


Hence ,  the  sequence 

(3.30)  (f0(p+m,‘)  €  L2(R0)  |  m  €  Z>  €  £  ®L2(RQ) 

mEZ 

and  the  operator  $  :  L2(B„)  *►  £  ®  L2(RQ),  defined  by 

0  ,p  mEZ 

(3.31)  $0>pf  -  (f0(p+m,-)  |  m  €  Z>  , 

is  an  isometry*  A  more  careful  application  of  the  Plancherel  theory 
shows  that  is  unitary.  Finally,  calculation  of  the  spectral  family 

{H0  p(li)  |  U  >  P2}  for  A0  gives 

(■(y-(p+m)2) 1/2 

(3.32)  H  (y)  f(X)  =  l  4»0  (X,p+m,q)  f  (p+m,q)  dq 

°’p  (p+m)2<UJo  ° 


In  particular,  making  p  ■*  <*>  gives  the  eigenfunction  expansion 


(3.33) 


f(X) 


L2 (B0)-lim 

M^oo 


0O (X  ,p+m,  q)  f0(p+m,q)  dq  . 


The  relationship  between  the  R-B  waves  for  A  and  the  reduced 
propagators  A^  will  now  be  discussed.  Note  first  that  if  lj;  is  an  R-B 
surface  wave  for  A  with  x-momentum  p  +  m  (-1/2  <  p  <  1/2,  m  e  Z)  and 
u)£(|p  +  £|  |  £€2}  then  €  D(A^0C)  and  for  y  >  h 
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(3.34)  <Hx,y)  ■  l  exp  {i(p+  4)x)  exp  {-((p+  l)1  -  w2) 1/2  y}  . 

|  p+JL  |  >to 

It  follows  that  <J)(x,y)  =  ij;(x,y)  |0  €  D(A  )  and  A  <J>  ■  U)2<£.  Thus  $  is  an 

**  P  P 

eigenfunction  of  A  .  To  formulate  the  converse,  note  that  every 

Hoc  Hoc 

<f>  €  L2  (ft)  has  a  unique  p-periodic  extension  iji  6  L2  (G)  .  It  is  easy 

to  verify  that  if 

(3.35)  =  ft  +  (2irn,0) 
then  for  each  m  €  Z  the  extension  ifj  is  given  by 

(3.36)  ip(x,y)  *  exp  {27rimp}  <p(x-2nm,y)  for  all  (x,y)  €  ft^  . 

This  defines  ip  in  L2°C(G)  because  G  differs  from  u„  by  a 

Lebesgue  null  set.  The  operator  0P:  L^oc(ft)  L2OC(G)  defined  by 
Hoc 

(3.36)  maps  L2  (ft)  one-to-one  onto  the  set  of  all  p-periodic  functions 
Hoc 

in  L2  (G) .  With  this  notation  it  is  not  difficult  to  show  that  if  <P  Is 
an  eigenfunction  of  then  ip  *  0^0  is  an  R-B  surface  wave  for  A  with 

reduced  x-momentum  p. 

The  relationship  between  R-B  diffracted  plane  waves  for  A  and 
generalized  eigenfunctions  of  A^  is  exemplified  by  (3.27).  More 
generally,  if  ^(X,p-hn,q)  is  an  R-B  diffracted  plane  wave  for  A  with 
-1/2  <  p  <  1/2,  m€  Z  then  0+(X,p+m,q)  »  \p±  (X,p+m,q)  |  satisfies 
<J>+(*  ,p+ra,q)  €  D(Ap0C),  (A  +  o>2  (p+m, q) )  <j>+  (X,p+m,q)  «  0  in  ft  and 

(3.37)  <t>±(X,p+m,q)  *  0Q  (X,p+m,q)  +  0  ’  (X,p+m,q)  ,  y  >  h, 


where  0^  (resp.,  0^)  has  a  Fourier  expansion  that  contains  only  outgoing 
(resp.,  incoming)  plane  waves  and  exponentially  damped  waves.  Functions 


4>+(X,p+m,q)  and  <f)_(X, p+m,q)  with  these  properties  will  be  called. 


respectively,  outgoing  and  incoming  diffracted  plane  waves  for  . 

They  are  unique  modulo  eigenfunctions  of  .  It  is  now  easy  to  verify 
that  if  <f>+(X,p+m,q)  (resp.,  $_(X,p+m,  q)  )  is  an  outgoing  (resp.,  incoming) 


diffracted  plane  wave  for  A^  then  ^+(X,p+m,q)  *  0P  $+(X,p+m,q)  (resp., 
l|/_(X,p+m,q)  -  0P  <(»_(X,p+m,q))  is  an  outgoing  (resp.,  incoming)  R-B 


diffracted  plane  wave  for  A  with  x-momentum  p  +  m.  These  relationships 
will  be  used  in  §6  to  construct  the  R-B  diffracted  plane  waves  for  A. 


§4.  Analytic  Continuation  of  the  Resolvent  of  A  . 


An  analytic  continuation  of  the  resolvent 
(4.1)  R(Ap>z)  -  (Ap  -  z)~l 

across  the  spectrum  0(Ap)  *  [p2,00)  is  constructed  in  this  section  by  an 
elegant  and  powerful  method  that  was  introduced  into  scattering  theory 
by  H.  D.  Alber  [3].  The  continuation  provides  the  basis  for  the  con¬ 
struction  in  §5  of  the  diffracted  plane  waves  <f>+(X,p,q)  for  Ap  and  the 
derivation  of  the  corresponding  eigenfunction  expansions. 

For  each  pair  of  extended  real  numbers  r,rf  satisfying  0  <  r 

<  r»  <  +oo 


B  .  -  (X  I 

1  -7r<x<TT,  r<v 

<  r'}  ,  B  =  B 

r,r 

‘ 

’  r  r," 

(4.2) 

ft  .  *  ft  n 

B  .  .  ft  *  ft 

r  ,r 

r  ,r 

r  r,« 

Moreover, 

!et  Pr  :  L2(Q0fr) 

-►  L2  (ft)  denote  the 

linear  operator  defined  by 

'  u(X)  , 

x  €  a 

o  ,r 

(4.3) 

Pr 

u(X)  *  < 

.  0  , 

X6ilr  . 

The  goal  of  §4  may  be  formulated  with  this  notation.  It  is  to  construct 
an  analytic  continuation  of 

(4.4)  z  -  R(Ap,z)  Pr  :  L2  (fl#  )  L2OC(ft) 

from  the  resolvent  set  p(Ap)  •  C  -  [p2,00)  across  cr(Ap)  -  [p2,«).  For 
this  purpose  p(Ap)  will  be  embedded  in  a  Rlemann  surface  M  . 
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The  definition  of  may  be  motivated  by  considering  the  linear 
space  of  functions 

(4.5)  Ep  z  r  ■  D(Ap°C)  n  (u  |  supp  (A  +  z)u  c  »  r  >  h  . 

Basic  properties  of  E  are  described  by 

r  P*z,r  7 

Lemma  4.1.  Every  u  €  E  satisfies 

-  P,z,r 

(4.6)  u  €  L*’£°C(nh) 

(4.7)  u(x,y)  -  l  u  (y)  ei^p+m^x  in  ft. 

n£Z 

where  the  series  converges  in  > 

(4.8)  um(y)  €  *  (h,«°)  . 

Moreover,  if  denotes  the  closure  of  in  R2, 

(4.9)  u  €  C°°(flr)  ,  and 

(4.10)  um(y)  »  c+  exp  {iy(z-  (p+m)2)1/2}  +  c~  exp  {-iy(z-  (p+m)2)1/2} 

for  y  >  h  where  c*  are  constants  and 
~  m 

(4.11)  Im  (z  -  (p  +  m)z)  >  0  . 

Properties  (4.6)  and  (4.9)  follow  from  elliptic  regularity  theory 

[1],  while  (4.7)  and  (4.8)  follow  from  classical  Fourier  theory.  The 

convergence  of  (4.7)  in  Lj’^oc(n^)  follows  from  the  fact  that  the  partial 

sums  of  the  Fourier  series  define  orthogonal  projections  in  L*((2  ,)  for 

2  r ,  r 

h  <  r  <  r’  <  ®.  (4.10)  follows  from  (4.9)  and  the  equation  Au  +  zu  ■  0 

in 
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Note  that  if  z  €  p(A  )  and  u  ■  R(A  ,z)P  f  with  f  €  L2 (0  )  then 

P  P  r  0  r 

u  €  L9 (ft)  n  E  and  hence  c  *  0  for  all  m  and  c+  *  0  when 

*  p ,  z ,  r  m  m 

Im  (z  -  (p  +  m)2)l//2  ■  0,  This  suggests  that  Mp  be  defined  as  the 
Riemann  surface  associated  with  the  family  of  holomorphic  functions  on 
C  -  [p2,00)  defined  by 

(4.12)  (z  (z  -  (p+m)2)1^2  |  Im  (z  -  (p+m)2)1^2  >  0  for  all  m  €  2}  . 


Mp  is  uniquely  determined  up  to  isomorphism  by  the  following  three 
properties  [3,  21]: 


(4.13)  Mp  is  connected  and  every  function  of  the  family  (4.12)  can 
be  continued  analytically  to  all  of  Mp. 

(4.14)  For  every  pair  of  points  of  Mp  that  lie  over  the  same  point 
of  C  there  are  at  least  two  functions  of  the  family  that  take 
different  values  at  these  points. 


(4.15) 


Mp  is  maximal  with  respect  to  these  two  properties. 


The  following  notation  will  be  used  in  connection  with  Mp.  £ 

will  denote  a  generic  point  of  M  and  tt  =  it  :  M  ■+  C  will  denote  the 

P  P  P 

canonical  projection  of  Mp  onto  C.  The  subscript  p  will  be  omitted  when 

there  is  no  danger  of  ambiguity.  The  analytic  continuation  of 

(z  -  (p  +  m)2)1/2  from  C  -  (p2,00)  to  M  will  be  denoted  by  wr.  (O  . 

p  p-ra 

Thus,  for  all  £  £ 


(A. 16) 


-  ±(t r(s)  -  (p  +  m)2)1/2 


M^  will  denote  that  component  of  Mp  over  C  -  [p2,00)  on  which 


Im  w (C)  >  0  for  all  m  G  Z.  Finally,  T»{(p  +  m)2  |  m  6  Z}  c  C  will 
p+m  p 

denote  the  set  of  branch  points  of  the  family  (4.12). 

The  properties  of  include  the  following.  has  infinitely 

many  sheets.  More  precisely,  for  each  disk  D(zQ,p)  C  C,  tt  l(D(zQ,p)) 

has  infinitely  many  components.  If  z0  *  (p+m)2  for  some  m  e  Z  then  the 

set  tt"1  (D(z0  ,p) )  contains  infinitely  many  branch  points.  Moreover,  for 

all  4  €  M  the  set  (m  I  Im  w  .(C)  <  0}  is  finite  [3].  Finally, 
p  p+m  - 

»  M  for  all  m  G  Z. 
p+m  p 

In  addition  to  M  the  set 

P 

(4.17)  M  -  |J  { (p»C)  |  5  e  Mp} 

-1/ 2<p<l/2 

will  be  needed  to  describe  the  dependence  of  the  continuation  of 

R(Ap,7T(C))Pr  on  p  and  C*  M  will  be  topologized  in  such  a  way  that  each 

function  (p,C)  w  .(C),  m  G  Z,  is  continuous  on  M.  To  this  end  let 
p+m 

(p0,C0)  e  M  and  define 

(4 . 18)  z0  =  it  (£„)  ,  D(z0,p)  =  {z  |  |z  -  z0|  <  p} 
and 

(4.19)  U(p0,£0>P>  *  ComP°nent  of  (D(zQ,p)) 

containing  C0  (C  M  ) 

Po 

To  define  a  neighborhood  basis  for  M  at  (pQ,C0)  three  cases  will  be 
distinguished. 

Case  1.  zQ  &  [p2,00).  If  p0  >  0  is  the  distance  from  z0  to 

[p2,oo)  then  for  p  <  p0  D(z0,p)  n  [pj,,00)  *  <$>  and  U(p0,C0>p)  contains  no 

branch  points  of  M  .  In  this  case 
*  o 
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(4.20) 


(sgn  Im  w  ^(O  |  m  €  Z}  ,  C,  e  U(p0,C0,p) 
P  0 


is  well  defined.  Moreover,  |p  -  pQ|  <6  implies  that  D(zQ,p)  n  [ p 2 , °°) 

■  <p  for  6  small  enough  and  hence  {sgn  Im  |  m  6  Z}  is  also  well 

defined  on  the  components  of  tt^1  (D (z0  ,p) ) .  In  this  case  one  may  define 
U(p,£0,p)  as  the  component  of  tt^1  (D(z0  ,p) )  for  which 

(4.21)  (sgn  Im  w (tt  *(z))  |  m  €  z}  =*  {sgn  Im  w  (tt~1(z))  |  m  6  Z) 
p-Hn  p  p0+m  p0 

for  z  6  D(z0,p).  A  corresponding  neighborhood  of  (p0,£  )  in  M  is 


defined  by 


(4.22) 


N(p0  ,S0  .P.'S)  *  U  Up, 5)  |  Z  €  U(p,C0,p)}  . 


|P-P0I<(S 


Case  2.  zQ  €  [p^,00)  -  T  In  this  case  if  p0  is  the  distance 

from  z0  to  the  set  T  then  for  p  <  p0  U(p  ,£0,p)  contains  no  branch 

Po  0  u 

points  of  M  and  (4.20)  is  well  defined  provided  tt  (?)  €  D+(z0,p) 

Pq  Po 

*  D(z0,p)  C\  {z  |  Im  z  >  0).  Moreover,  |  p  -  pQ  |  <  6  implies  that  D(zQ,p) 

contains  no  points  of  T  ,  for  6  small  enough,  and  hence 

{sgn  Im  w_.  (C)  I  in  €  Z}  is  also  well  defined  if  tt  (5)  e  D,(zn,p).  In 

pH’Ul  P  T  " 

this  case  one  defines  U(p,£0,p)  as  the  component  of  tt”1  (D(zQ  ,p) )  for 
which  (4.21)  holds  for  z  €  D+(z0,p)).  A  corresponding  neighborhood  is 


again  defined  by  (4.22). 


Case  3. 


(p0  +  m0)2  for  some  mQ  €  Z.  If  p.  >  0  is  the 


distance  from  z0  to  the  set  T  -  {(pQ  +  mQ)2}  then  for  p  <  p0  the  set 

Po 

u(Po*Cn>p)  contains  only  one  branch  point;  namely,  that  for  w  _  (c)  • 

Po+*o 

Hence  {sgn  Im  w  (<;)  I  m  €  Z  -  {mrt}}  is  well  defined  for 
p0+m  0 

5  €  U  (p0  ,?0  ,p)  and  tt  (O  €  D  (zQ,p).  Moreover,  |  p  -  p0 1  <6  implies 

*  0 


that  D(z0,p)  contains  (p  +  mQ)2  and  no  other  points  of  the  set  and 
hence  {sgn  Im  |  m  G  Z  -  {mQ}}  is  well  defined  on  the  components 

of  (D+(zQ  ,p))  .  In  this  case  one  may  define  U(p,<;o,p)  as  the 
component  of  (D(zQ  ,p) )  for  which 


(4.23) 


(sgn  Im  wp+m(7rp1  (z))  I  ro  e  Z  -  {m0}} 


■  {sgn  Im  wp  ^  ( TTp ^  (  2  )  )  |  m  6  Z  -  {mQ}} 


for  all  z  G  D+(z0,p).  A  corresponding  neighborhood  is  again  defined  by 


(4.22), 


The  topology  of  M  is  defined  to  be  the  one  generated  by  the 


neighborhood  bases  defined  above  and  one  has 

Theorem  4.2.  Each  of  the  functions  on  M  defined  by 


(4.24) 


( Pf Z)  +  Wp-Hn(^  >  tn  6  Z 


is  continuous  on  M.  Moreover,  the  family  of  functions 

*  Vm<«  I  m  G  z}  is  equicontinuous  in  M. 

The  theorem  that  iz  w^(£)  I  m  6  Z}  is  equicontinuous  on  M 

p+m  p 

for  fixed  p  was  proved  by  Alber  [3].  Theorem  4.2  plays  a  key  role  in 

proving  the  continuity  in  (p,q)  of  the  Rayleigh-Bloch  waves  in  §6. 

The  Frechet  Space  F  „  .  To  describe  the  subset  of  E  that 

_  _  _  ■  PtC>r  P,z,r 

contains  the  analytic  continuation  of  R(A^,z)Prf  to  consider  the  set 

of  functions  uG  E  whose  Fourier  representations  (4.7),  (4.10)  in 

P  >  z  >  r 

&r  satisfy 


(4.25) 


For  each  m  s  Z,  either  c  *  0  or-c  *  0,  and 

m  m  -  - 


(4.26) 


0  for  all  but  a  finite  number  of  m  €  Z  . 
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Note  that  these  conditions  express  the  "radiation  conditions" 


(4.27) 


[0^  ±  i(z  -  (p  +  m) 2 ) 1/2  ]  um(y)  *  0  ,  y  >  r  , 


where  for  each  m  either  "+"  or  is  chosen  and  is  chosen  for  all 

but  a  finite  number  of  m  G  Z.  It  is  clear  that  each  such  u  G  E  is 

pfz,r 

associated  with  a  unique  point  £  G  such  that  tt  (£)  *  z  and  the  Fourier 
expansion  (4.7),  (4.10)  of  u  has  the  form 


(4.28) 


u(x,y)  =  l  c  exp  {i(p+m)x  +  iyw  (?)}  ,  y  >  r 
mEZ  m  p+m  - 


For  each  (p,£)  e  M  and  each  r  >  h  the  set  of  all  such  solutions  will  be 
denoted  by 


(4.29)  F  =  D(A^0C)  n  {u  |  supp  (A+tt  (£))u  C  Q.  and  (4.28)  holds} 

P>  4 >  r  P  P  o ,  r 

Note  that  F  C  D(A^°C)  C  L*’^°c(A,fl)  and  recall  that  D(A^OC)  is 

p,4»r  P  2  P 

closed  in  L* ’^OC (A,ft) .  This  implies 

Theorem  4.3.  F  is  closed  in  D(A^°C)  in  the  topology  of 

P  >  4 » r  P 

L’AoC(A,n)  and  hence  is  a  Frechet  space. 

This  is  immediate  because  the  defining  properties  of  ^  ^ , 

namely  supp  (A  +  it  (£))u  C  fl  and  [D  -  i  w  (?)]  u  «  0  in  y  >  r, 
p  o  ,r  y  p-rm  m  — 

are  preserved  under  convergence  in  L2*^0C(A,ft). 

The  following  condensed  notation  will  be  used  in  discussing 

F  and  related  operators: 

P  9  4  > r 

(4.30)  (u,v)  ,  *  (u,v). 


r  ,r 


Mnr,r'>  ' 


(4.30) 
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(4.31) 

<u-v\;r,r'  *  (u‘v)lJ «r  r ,)  • 

(4.32) 

<U’V)l;4!r,r'  •  <u'v)Lj(A,!Jt  r,)  . 

Now  let  P  _  : 

P>5»r 

F  -►  L2  (Si  )  denote  the  natural  projection  defined  by 

P)Si^  0  >  ^ 

(4.33) 

P  _  u  =  u  I  n  for  all  u  e  F  _ 

p,C,r  p,C,r 

An  important  property  of  F^  ^  is  expressed  by  the  following  generali¬ 
zation  of  a  theorem  of  Alber  [3,  p.  264], 

Theorem  4,4,  For  every  compact  set  K  C  M  and  for  every  r*  >  r 
there  exists  a  constant  C  =  C(K,r,rf)  such  that 


(4.34) 

for  all  u  €  (^J 
isomorphism  of  F 


lull 


a  f  <  c  a  p  „  u# ,  A 

i ; A; o,rT  -  p>C»r  i ;A; o ,i 


/  •  \/-*Tr  F  „  .  In  particular,  P  is  a  topological 

(p,C)6K  p,£,r  v  *  p,£,r  *  & 

„  onto  P  F  ,  topologized  by  the 

p,£,r  p,£,r 


Lj(A,Q  )-norm. 


o,r 


The  Operators  Ap^r  :  L2  L2  *  Following  Alber1  s 
program,  the  construction  of  the  analytic  continuation  of  R(Ap,z)Pr  to 
M  will  be  based  on  the  family  of  linear  operators  A 


*  _  in  L- (Q  ) , 

p,C,r  2  o ,r' ’ 


defined  for  all  (p,C)  €  M  by 


(4.35) 

(4.36) 


D(A  )  -  P  F  , 


A  u  ■  -Au  . 

P»C,r 


The  properties  of  Ap  ^  r  that  are  fundamental  for  the  analytic  continu¬ 
ation  of  R(Ap,z)  are  described  by  the  following  theorems. 
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Theorem  4.5.  For  every  (p,£)  €  M  and  every  r  >  h  the  operator 
A  is  m-sectorial  in  the  sense  of  Kato  [17,  p.  279]. 

P  *  S  9  r 

Theorem  4.6.  For  all  grating  domains  of  the  class  defined  in  §1, 

the  family  of  operators  {A  |  (p>C)  e  M}  is  continuous  in  the  sense 

of  generalized  convergence  (Kato  [17,  p.  206]).  Moreover,  for  each 

fixed  p  €  (-1/2, 1/2]  the  family  {A  ]  C  €  M  }  is  holomorphic  in  the 

generalized  sense  (Kato  [17,  p.  366]). 

Theorem  4.7.  For  every  (p,£)  6  M,  every  r  >  h  and  every 

z  €  p(A  )  the  resolvent  R(A  ,z)  =  (A  -  z)  1  is  a  compact 

P  >  4 »  r  P »  4  > r  P » 4  >  r 

operator  in  L2  (Qq  ^)  and  hence  a(Ap  ^  r)  is  discrete. 

Theorem  4.5  generalizes  Alber  [3,  Th.  5.5].  As  in  [3]  it  may  be 

proved  by  associating  A  with  a  densely  defined,  closed,  sectorial 

P  >  4 » r 

sesquilinear  form  in  L2 (ftQ  ^)  and  using  Kato's  first  representation 
theorem  [17,  p.  322].  The  second  statement  of  Theorem  4.6  generalizes 
Alber  [3,  Th.  5.5b].  The  hypothesis  G  6  S  of  §1  is  needed  to  prove 
Theorem  4,6.  Theorem  4.7,  which  generalizes  Alber  [3,  Th.  5.5a],  is  a 
consequence  of  the  local  compactness  property  of  G  in  the  case  of  the 
Neumann  boundary  conditions.  Complete  proofs  of  Theorems  4.5,  4.6  and 
4.7  are  given  in  [38].  The  following  consequences  of  these  theorems  are 
needed  for  the  spectral  analyses  of  Ap  and  A  in  §§5-6. 

Theorem  4.8.  For  all  £  €  M*  one  has  77^(0  G  p(A^  £  r)  an<* 


(4.37) 


’t<V5,r'V«)  '  V?,r  »‘VV5))  Pr 


This  result  may  be  verified  by  direct  calculation. 
Theorem  4.9.  For  every  p  €  (-1/2, 1/2]  the  set 


(4.38) 


it  (;)  e  a  (A  )}  c 
P  P»?.r  P 
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has  no  accumulation  points  in  and  is  independent  of  r  >  h. 

This  result,  which  generalizes  [3,  Th.  5.5c],  is  a  consequence 
of  Theorem  4.7.  For  brevity  the  resolvent  of  (4.37)  will  be  denoted  by 


(4.39) 


R  r  „  -  R(An  »ir  (C))  €  B(L 2  (ft  )) 
p,5,r  p,5,r  p  o  ,r 


Here  B(X)  denotes  the  bounded  operators  on  X. 

Corollary  4.10.  For  each  p  €  (-1/2, 1/2]  and  r  >  h  the  mapping  * 


(4.40) 


C  +  R  €  B(L2(ft  )) 
p,C,r  o,r/y 


is  finitely  meromorphic  on  M  with  pole  set  Z  . 

P  P 

This  result  is  based  on  a  theorem  of  S.  Steinberg  [28];  cf. 

[3,  Th.  5.5e].  Theorem  4.4  and  4.8  provide  the  analytic  continuation  of 
R(Ap,z)  Pr  in  the  following  form. 

Corollary  4.11.  The  analytic  continuation  to  of 

(4.41)  c  ♦  R(Ap,irp(C))Pr  €  B(L2(n0>r)  ,  L2,JloC(A,0))  ,  C  €  it  , 


is  given  by 


(4.42) 


'JU.r  V<,r  *  *«•*<«„, r>>  ‘I'10'11'®''  «  e  «, 


where  B(X,Y)  denotes  the  bounded  linear  operators  from  X  to  Y. 

Corollary  4.12.  For  all  grating  domains  of  the  class  defined  in 
Si,  the  point  spectrum  a0(Ap)  is  discrete. 

This  result  follows  from  Theorem  4.9  and  Corollary  4.10. 
Corollary  4.13.  For  all  grating  domains  of  the  class  defined  in 


§1  one  has 


ii  ii  ii.  ii 
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(4.43) 


n  (AT-  n  l  )  c  on  (A  )  u  T 

p  p  P  o  p  p 


where  is  the  closure  of  M+  in  M  . 

P  P  P 

If  CJ0(A  )  *  <|>  then  (4.43)  means  that  the  poles  of  R  that  lie 

above  the  spectrum  a(A^)  *  [p2»°°)  must  lie  above  the  branch  point  set 

T  .  The  fact  that  such  poles  may  or  may  not  occur  is  illustrated  by  the 
D  N 

two  operators  AQ  and  AQ  corresponding  to  the  degenerate  grating.  For 
A^,  separation  of  variables  leads  to  a  construction  of  the  Green's 
function  (*  kernel  of  the  resolvent  R(A^,z))  which  can  be  written 


G^(X,X\p,z) 

(4.44) 

-4  l  ei(P^)(x-x’)(2_(p^)2)-I/2sin  Cz_Cp+m)2),/2  y<  ei<-(P*0*)1'*  y> 
m€Z 


where  y<  *  Min  (y,yf),  *  Max  (y,yf).  The  analogous  calculation  for 

A0  gives 

Gj(X,X’,p,z) 

(4.45) 

i  r  i(p+m)(x-xr)  /  ,  ,  v2\-i/2  /  /  \  2  \  1/2  i  ( z-  (p+ra) 2 ) 

■xr  l  e  (z-  (p+m)  )  cos  (z-(p-Kn) z) y  e  r  '  7> 

1  roEZ 

In  the  first  case  R(A^,z)  has  no  poles  for  real  z  ■  A  ±  iO  €  [p2,«).  In 
N 

the  second  case  R(A0,z)  has  a  simple  pole  at  each  of  the  points 

z  -  X  ±  iO  €  T  . 

P 

The  following  two  theorems  are  implied  by  Theorems  4.4  and  4.6. 


m 
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Theorem  4.14.  Let 


2  -  U  Up,;)  |  ;e  £  }  ; 


-p.  U  {(p,;)  I  TTp(;)€a(Ap>i;>r)} 

-1/ 2<p<l/2  -l/2<p<l/2 


(4.46) 


Then  M-I  is  open  in  M  and 


(4.47) 


<»•«  *  %,C,r  S  B<WW.,r» 


is  continuous  on  M-Z. 

Theorem  4 >15.  The  mapping 


(4.48) 


(P,0  -*■  P'1  R 


6  B(L2<no,r)>L2,£°C(A’n)) 


is  continuous  on  M-E. 

A  direct  consequence  of  Theorem  4.15  that  is  needed  below  is 

% 

Corollary  4.16.  Let  K  be  any  compact  subset  of  M-L  and  let 
r’  >  r  >  h.  Then  there  exists  a  constant  C  =  C(K,r,r')  such  that 


(4.49) 


IP. 


-l 


fll 


p,£,r  p,;,r  i;A;o,r  - 


.,  <  C  If! 


for  all  (p,£)  e  K  and  all  f  e  L2(fiQ  f) . 

A  Limiting  Absorption  Theorem.  In  the  remainder  of  this  report 

the  point  £  will  be  restricted  to  M+,  the  closure  of  M"1"  in  M  .  To 

P  P  P 

simplify  the  notation  points  ;  €  Mp  will  be  identified  with  their  images 
TTp(C)  ■  2  e  C  -  [p2,®)  and  the  points  of  3Mp  will  be  denoted  by  X  ±  iO, 
where  X  e  f p 2 , °°)  -  With  this  notation  the  operators 


C±io,r  Rp.U,a.rs  >0',<»..r)'Li,*OCW-a» 


(4.50) 
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are  defined  and  continuous  for  all  X  ±  io  €  Af1’  -  Z  •  Note  that  by 

P  P 


Corollary  4.13 


(4.51) 


o(A)  -  a.  (A  )  -  T  c  i(3<  -  E  )  . 


P  P 


Now  let  f  €  L-(Q  )  and  define 
2  o  *r 


-i 


(4.52)  U.C.P.A)  -  PpiXll0,t  «p,Ul.,r  f  €  Fp>xtl,>r 

Then,  in  particular, 


(4.53) 

(4.54) 


ioc. 


u±(*,p,X)  €  D(Ap  )  ,  and 


A  u+  +  X  u+  *  f  in  0  . 


Moreover,  tt  (X  ±  iO)  -  X  for  all  X  €  a(A  )  and 
P  p 


wpfm(X  ±  iO)  -  ±(X  -  (p  +  m)2)l/2  if  X  >  (p  +  m)2  , 


(4.55) 


Wp-Hn^  ~  *  i((p  +  m)2  -  X)"2  if  X  <  (p  +  m)2  . 


Hence,  the  Fourier  series  (4.28)  of  u+  have  the  form 


u±(x,y,p,X)  -  l  c*  ei(p4m)x  *) 


2  \  1/2 


(p+m) 2<X 


(4.56) 


+  c  ±  ei(p4in)x  -y((p+m)2-X) 
4  m 


(p+m)2>X 


Thus  u+  and  u_  are  the  outgoing  and  incoming  solutions,  respectively,  of 
the  boundary  value  problem  (4.53),  (4.54).  Moreover,  they  are  uniquely 
determined  by  these  conditions,  by  Theorem  2.1,  provided 


-  ■  — 


t 
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(4.57)  X€  a(Ap)  -  a0(Ap)  -  Tp  . 

The  final  result  of  this  section  is  a  uniform  bound  for  the 
functions 


(4.58) 


Pn\-Mrr  Rn  *  6 

P ) A— io ) r  p i A- io j r 


which  may  be  formulated  as  follows. 

Corollary  4.17.  Let  I  *  [a,b]  satisfy 


(4.59)  ICO(Ap)-Ot(Ap)-Tp 

and  let  pf  aQ,  r  and  r'  satisfy  -1/2  <  p  <  1/2,  a0  >  0  and  r'  >  r  >  h. 
Then  there  exists  a  constant  C  ■  C(I,p,a0 , r ,r ' )  such  that 


(4.60) 


IP-*  .  R  .  .  .  fl  .  ,  <  C  Ifl 

p,X±ia,r  p,A±xa>r  i;A;o,r  -  o,r 


for  all  A  €  I,  0<a<a0  and  all  f  e  L2(H  ). 

—  —  0  >  r 

This  result  is  a  direct  consequence  of  Corollary  4.16. 


J 


a 


§5,  The  Eigenfunction  Expansions  for  A  . 

_ E 

This  section  presents  a  construction,  based  on  the  limiting 
absorption  theorem  of  §4,  of  the  diffracted  plane  wave  eigenfunctions 
4>+(X,p+m,q)  and  a  derivation  of  the  corresponding  eigenfunction 
expansions  for  Ap.  For  brevity  the  derivation  is  restricted  to  the 
cases  for  which  0o(Ap)  -  <j>.  The  modifications  that  are  needed  when 
ap(Ap)  +  4>  are  indicated  at  the  end  of  the  section. 

Throughout  this  section  p  €  (-1/2, 1/2]  is  fixed,  m  G  Z  and  q  >  0. 
<t>0  (X,p4m,q)  denotes  the  generalized  eigenfunction  for  Aq  ;  that  is,  one 
of  the  functions  (3.25),  (3.26).  The  corresponding  outgoing  and  incoming 
diffracted  plane  waves  for  Ap  are  characterized  by  the  properties 

(5.1)  <M*,P-Hn,q)  €  D(ApOC)  , 

(5.2)  (A  +  0J2(p+m,q))  <f>+(X,p+m, q)  =  0  in  ft  , 

(5.3)  <j>+(X,p*+tn,q)  -  (X,p-kn, q)  +  $j(X»p4m,q)  ,  y  >  h  , 

where  (resp.,  <^)  is  an  outgoing  (resp.,  incoming)  diffracted  plane 
wave  in  These  properties  imply  the  symmetry  relation 

(5.4)  4>.(X,p+m,q)  -  (J)+(X,-p-m,q)  . 

Hence  it  will  be  sufficient  to  construct  the  functions  4>+(X,p+m,q)  . 

To  construct  <J>+  fix  an  r  >  h  and  introduce  a  function 
j  €  C°°[0,«)  such  that  j  f  (y)  >0,  0  <  j  (y)  <  1,  j  (y)  =  0  for 
0  <  y  <  (h  +  r)/2  and  j (y)  =  1  for  y  >  r.  Next  define  the  function 
<fc+(X,p+m,q)  for  all  X  €  Q  by 
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(5.5)  4>+(X,p-hn,q)  *  j  (y)  <p0  (X,p-Hn,q)  +  <J>^(X,p+m,q)  ,  X  e  ft  . 

Then  (5.1),  (5.2),  (5.3)  imply  that  <f>_|_  is  characterized  by  the  properties 

(5.6) 


<f>|(*,p+m,q)  €  D(Ap°C) 


(5.7) 


(5.8) 


(A  +  u>2(p-hn,q))  <|>^(X,p+m,q)  *  -M(X,p+m,q)  in  ft  , 
^(X,p+m,q)  is  an  outgoing  diffracted  plane  wave. 


The  function  M  in  (5.7)  is  defined  for  all  X  €  R^,  p+meR  and  q  >  0  by 


(5.9) 


M(X,p+m,q)  =  (A  +  co2(p-hn,q))  j  (y)  <f>0 (X,p+tn,q) 

=  j"(y)  <J>0(X,p+m,q)  +  2  j’(y)  D2  $0 (X,p+m,q) 


and  has  the  properties 

(5.10)  M  6  C°°(R2  x  R  x  Rfl)  , 

(5.11)  M(x+2n,y,p+m,q)  =  exp  {2rrip}  M(x,y,p+m,q)  , 

(5.12)  supp  M(*,p+m,q)  C  {X  |  (h  +  r)/2  <  y  <  r}  . 

It  follows  that  M(* ,p+m,q) |n  £  L2(2  )  and  hence  (5.6),  (5.7),  (5.8) 

wo,r  ,r 

can  be  integrated  by  means  of  the  analytic  continuation  of  the  resolvent 
of  Ap  defined  by  (4*50).  More  generally 

(5.13)  <J>’  (•  ,p+m,q,z)  -  P"1  R  M(*,p+m,q)  €  D(a5°C) 

P»z,rp,z,r  p 

and  z  +  (pr  ( •  ,p+m, q, z)  6  L2*^°C(A,ft)  is  continuous  for  all  q  >  0  and 
z  €  Mp  -  £p.  Hence,  4>f  (•  ,p-Kn,p, A+iO)  €  D(Ap0C)  satisfies  (A  +  ■  0 

in  ft  and  the  outgoing  radiation  condition  (4.56)  for  all  X  €  [p2»°°)-Tp. 


t 

i 


i 

! 


O 
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In  particular,  the  solution  of  (5.6),  (5.7),  (5.8)  is  defined  by 

(5.14)  0^(*,p-hn,q)  «  <)> 1  (•  ,p+m,q,w2  (p+m, q)  +  iO) 

for  all  q  €  R  -  E  where 
M  0  m,p 

(5.15)  Em  p  -  (q  >  0  |  u)2(p+m,q)  E  Tp}  . 

Note  that  E  is  a  countable  subset  of  R0  *  (0,°°)  with  no  finite  limit 
m,p 

points. 

The  diffracted  plane  wave  <t>+(X,p+m,q)  is  defined  by  (5.5), 
(5.13)  and  (5.14)  and  one  has 

Theorem  5.1.  Let  G  be  a  grating  domain  of  the  class  defined  in 
§1  and  let  a0  (Ap)  *  <j>.  Then  there  exist  unique  diffracted  plane  wave 
eigenfunctions  4>+  (X,p+m, q)  for  each  p  E  (-1/2, 1/2],  m€  Z  and 

j  £oc 

q  €  R0  -  E  .  Moreover,  q  -►  <J>+(*  ,p+m,q)  E  L2  9  c(A,ft)  is  continuous 
m,p  — 

for  q  E  Rrt  -  E 

0  m,q 

The  uniqueness  follows  from  Theorem  2.1  and  aQ  (Ap)  «  <p.  The 
continuity  is  a  consequence  of  Theorem  4.15. 

The  functions 

J?  Ap 

(5.16)  <J>(X,p+m,q,  z)  »  j  (y)  <pQ  (X,p+m,q)  +  <f> ’  (X,p+m, q,  z)  E  D(Ap  )  , 

which  are  defined  for  p  E  (-1/2, 1/2],  m  E  Z,  q  >  0  and  z  e  M+  -  I  will 

P  P 

be  used  in  deriving  the  eigenfunction  expansions  for  <j>+  and  <J>_.  They 
will  be  called  approximate  eigenfunctions  of  Ap  because 

(5.17)  (A  +  z)  <j>(X,p+m,q,z)  •  (z  -  u)2(p+m,q))  j  (y)  4>0(X,p+m,q) 


and 


Then  IJI  ■  1  and  (5.21)  implies 


(5.23) 


IIR(Ap,z)ffi  2  »  II J  R(Ap ,  z)  f#  2  +  0(1)  ,  In  z  +  0  , 


uniformly  for  Re  z  S  I.  Next,  Parseval's  relation  (3.29)  for  A #  and 

(5.23)  imply 


(5.24)  lR(An,z)f»2  -  J  I  <J  R(Ap,z)f) ~  (p+m,-)#2  +  0(1),  Im  z  -*•  0, 


uniformly  for  Re  z  €  I.  To  relate  this  to  the  eigenfunctions  <p+  define 


(5-25) 


and  note 


f(p+m,q,z)  =  <J>(X,p+m,q,z)  f(X)  dX,  f  €  L2  (ft) , 

'ft 


Lemma  5.2.  For  all  f  €  L,  (ft)  one  has 


(5.26) 


f  (p+m,q,  z)  »  (aj2(p+m,q)  -  z)  (J  R(A  ,z)f)~  (p+m,q)  . 


A  heuristic  proof  of  (5.26)  is  contained  in  the  following  formal 
calculations,  based  on  (5.17). 


(5.27) 


f(p+m,q,z)  »  R(A  ,  z)  (A  -  z)  <|>(X,p-Hn,q,  z)  f(X)  dx 
Jo  P  P 


j  (to2  (p+m,q)  -  z)  j  (y)  <pg  (X,p+m,q)  R(A  ,  z)  f(X)  dX 
Jo  P 


(ur(p+m,q)  -  z)  <J>0(X,p+m,q)  j  (y)  R(A.z)  f(X)  dx 


(uj2 (p-Hn, q)  -  z)  (J  R(A  ,z)f)T  (p+m,q)  . 


The  calculation  is  not  rigorous  because  the  presence  of  the  term  jif0  in 
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(5.16)  implies  that  <\>(  •  ,p+m,  q,  z)  ^  DCA^),  A  rigorous  but  longer  proof 
may  be  given  by  the  technique  of  [34,  p.  94]. 

Combining  (5.24)  and  (5.26)  gives 


(5.28) 


nR(Ap,z)fll2 


l 

mEZ 


f (p+m, * , z) 
0)2(p+m,  •  )-z 


2 

+  0(1)  ,  Im  z  *►  0, 


uniformly  for  Re  z  G  I.  Hence,  putting  z  =  A  ±  i a,  multiplying  by  q/tt 
and  integrating  over  A  E  I  gives 


a 

7T 


'i 


HR(Ap,A±ia)fll2 


I 

mEZ 


f  (p-Hn,q,A±ia)  1  2 
(  A— qj 2  ( p+m ,  q)  )  2+C2 


dq  dA  +  0(a) 


(5.29) 


•  00 


=  ^ 
m£Z  '  o 


£  f  |  f  (p-hn,q,A±ia)  |  2  . 

,TT  Jx  (A-(jl)2  (p+m,  q)  )  2+a2  aA 


dq+  0(a) 


by  Fubini's  theorem.  The  determination  of  11^(1)  will  be  completed  by 
calculating  the  limit  for  a  +  0  of  the  last  equation.  Note  that  the 
continuity  of  the  approximate  eigenfunctions  (5.16)  for  q  >  0, 
z  E  -  Xp  (cf.  (5.13))  implies  that  f (p+m, q, A±ia)  is  continuous  for 
q  >  0,  A  E  [p2,00)  -  T  ,  a  >  0.  Thus  if  one  defines 

(5.30)  f  (p+m,q)  =  f  (p+m,q,0)2(p+m,q)  +  iO)  ,  q€  R  -  E 

-  u  m,  p 

then  for  all  f  €  L.2°m(£2) 


(5.31) 


f±(p+m,q) 


<j)+(X,p+m,q)  f(X)  dx 
£2  " 


and 


(5.32) 


f±(p+m. •)  e  C(R0 


E 

m,p 


)  • 
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The  calculation  of  the  limiting  form  of  (5.29)  will  be  based  on  the 
following  two  lemmas. 

Lemma  5.3.  For  every  f  £  and  every  closed  interval 

I  C  [p2,00)  “  Tp  one  has 


(5.33) 


lim  “ 
0+o  + 


[  [f  (pHa,q,A±ig)  |2  A 
Jj  (A-U)2  (p+tn,q) )  2+a2 


XI(w2(p-hn,q))  |f_(p+m,q)  |2 


for  all  q  e  Rn  -  E  where  Xt(^)  is  the  characteristic  function  of  I, 
m,p  l 

normalized  so  that  Xj(a)  =  Xj-Cb)  =  1/2 . 

Lemma  5.3  follows  from  the  continuity  of  f  (p-Hn,  q,  X±ia)  and  well- 
known  properties  of  the  Poisson  kernels;  cf.  [34,  p.  101]. 

Lemma  5.4.  For  every  f  G  L^^ft)  ,  every  closed  interval 
1  C  [p2,00)  -  Tp  and  every  a0  >  0  there  exists  a  constant  C  =  C(f,I,afl) 
such  that 


(5.34) 


l  [  |f (p+m, q, Alia) | 2  dq  <  C 
mEZ  '  o 


for  all  p  €  (-1/2, 1/2],  X  6  I  and  a  e  [ 0 , a0 ] . 

This  result  is  the  analogue  of  [34,  Lemma'  6.8,  p.  103].  A  full 
proof,  based  on  Corollary  4.17,  is  given  in  [38]. 

The  limit  of  equation  (5.29)  for  a  -*■  0  may  now  be  calculated. 
Lemma  5.3  gives  the  limits  of  the  inner  integrals  in  (5.29).  Term-wise 
passage  to  the  limit  can  be  justified  by  Lemma  5.4  and  Lebesgue's 
dominated  convergence  theorem;  see  [34]  for  details.  The  result  is, 
by  (5.19), 


»np(I)fll2 


.00 

l  XrO'^Cp+m.q))  |f+(p-Hn,q) 

tn=Z  J 


dq 


(5.35) 


0 
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for  all  f  €  L2°m(ft)  and  I  c  [p2 ,°°)  -  Tp  where  f+(p+m,q)  is  given  by 
(5*31) . 


The  Eigenfunction  Expansions  for  A^.  The  eigenfunction  expan¬ 
sions  for  Ap  based  on  ^(Xjp+m, q)  and  4>_(X,p+m,q)  can  be  derived  from 
(5.35)  and  the  spectral  theorem  by  standard  methods;  cf.  [34,  p.  109ff] 
Only  the  results  are  given  here.  Details  may  be  found  in  [34]  and  [38] 


To  begin  note  that  since  a0 (A^)  *  $  the  restriction 


I  C  [ p z , °°)  -  Tp  can  be  dropped;  (5.35)  is  valid  for  f  €  L2°m(^)  and  all 


I  C  [p  ,°°) .  Making  I  +  [p  ,°°)  then  gives  the  Parseval  relation 


(5'36)  "f"L2(fi)  =  Jz  »'f±(p+m«,)^2(R0) 

for  all  f  €  L2°m(0).  Together  with  (5.32)  this  implies  that  for 
f  e  L^om(S2)  , 

(5.37)  f±(p*n,-)  6  C(R0  -  E  )  n  L2(R„). 

A  standard  density  argument  then  implies 

Theorem  5.5.  For  all  f  e  1.2(0)  the  limits 

(5.38)  f+(p+m,q)  =  L2(R0)-lim  [  <f>+(X,p+m,q)  f(X)  dX 

exist  and  (5.35),  (5.36)  are  valid  for  all  f  €  L2  (ft)  . 

An  eigenfunction  representation  of  the  spectral  family  can  now 
be  obtained  from  (5.35)  by  the  usual  polarization  and  factorization 
arguments.  In  this  way  one  obtains 

Theorem  5.6.  For  all  f  e  L2(ft)  one  has 

r(U“(P+ni)2)  1/2 

(5.39)  n  (y)  f(X)  -  l  <t>+(X,f*m,q)  f+(p+m,q)  dq 

w  (p+m)2<y  J  o 
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and  hence 


(5.40) 


f(X)  «  L2 (fi)-lim 
M-«o 


fM 

,  I  ♦* 

|  m J  <M  J  o 


(X,p+m,q)  f+(p+m,q)  dq  . 


Finally,  define  linear  operators 


(5.41) 

by 


:  L2(R)  -  l  ®  L2(R0) 

~,P  nez 


(5.42) 


$+  f  *  ( f ^ (p-Hn,  • )  |  m  €  Z) 
- » p  - 


Then  $  and  $  are  spectral  mappings  for  A  in  the  sense  of 
“**»P  "  9  P  P 

Theorem  5.7.  For  every  bounded,  Lebesgue-measurable  function 
'F(A)  defined  on  p2  <  X  <  00  one  has 


(5.43) 


(♦*, p  '*rCAp)f)m  =  'i'(0)2(p^,-))($±>pf)m  ,  m€  Z 


where  ^(A^)  is  defined  by  the  spectral  theorem. 

Finally,  the  orthogonality  and  completeness  of  the  generalized 

eigenfunctions  (f>+  is  expressed  by 

Theorem  5.8.  The  operators  and  <t>_  are  unitary. 

It  Is  clear  from  Parsevalfs  relation  (5.36)  that  are 

±»P 

isometries  which  proves  the  completeness  relation 


(5.44) 


* 

<t>  $  »  1 

+  n  *  4-  rt  A 

-  »r  -  *  P 


The  surjectivity  of  which  is  equivalent  to  the  orthogonality 

x ,  p 

relation 


$  $  «  1 
±>P  *±,P 


(5.45) 
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is  not  a  consequence  of  the  spectral  theorem.  A  proof  of  (5.45)  by  the 
method  introduced  in  [34,  p.  112ff]  is  given  in  [38], 

Operators  that  have  Point  Spectrum.  It  was  shown  in  §4  that, 
in  general,  aQ(Ap)  is  discrete.  Let  5C0  be  the  subspace  of  L2(ft)  spanned 
by  the  eigenvectors  of  A^  and  let  dim  JC0  *  N(p)-  1  <  oo#  Let 
(^j(p)  |  1  f  j  <  N(p)}  be  the  eigenvalues,  repeated  according  to  their 
multiplicity  and  enumerated  so  that  Xj (p)  <  \.+1 (p) .  Let 
{^j (X,p)  |  1  <  3  <  N(p)}  be  a  corresponding  orthonormal  set  of 
eigenfunctions . 

Proceeding  as  before  it  is  found  that  the  diffracted  plane  waves 
([)+ (X,p+m,  q)  can  be  constructed  and  Theorem  5.1  holds  with 

(5.46)  ^m,p  =  >  0  I  w2(p+“.q)  €  Tp  U  a0(Ap)} 

which  is  still  a  countable  set  with  no  finite  limit  points.  Similarly, 
the  spectral  family  (11^  (y ) }  still  satisfies  (5.35)  for  f  E  L2°m(ft)  if 
I  C  [p2,00)  -  Tp  -  a0(Ap).  It  follows  that  II  (y)  differs  from  (5.39)  only 
by  the  projection 


(5.47) 


l  <t>,  (X,p)  f,(p)  ,  f,(p) 
Xj (p)<y  J  J  J 


(4>j  (*  .p)  *f >La  • 


and  Parseval's  relation  and  the  eigenfunction  expansion  become 


N-l  „ 

(5.48)  0  f  1 2  *  l  |f.(p)|2  +  l  I  f+(p-hn,  •)!  2  ,  f  €  L2  (Q)  . 

j-1  J  meZ 


and 


N-l 

(5.49)  f(X)  -  l  4i.(X,p)  f.(p)  +  l  <J>+(X,p-Hn,q)  f  (p+m,q)  dq  , 
j-1  2  J  n£Z  J0  ‘ 
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convergent  in  L2(ft).  The  form  of  the  spectral  family  implies  that  Ap 
has  no  singular  continuous  spectrum  :  L2  (ft)  *  3C0  ®  3C  ,  where  is  the 

aw  a  v 

subspace  of  absolute  continuity  for  [17,  Ch.  X].  Finally,  Theorem 

5.8  must  be  modified  to  state  that  and  <t>  are  partial  isometries 

+»P  “»P 

with  initial  set  and  final  set  I  &  L2(RQ): 


(5.50) 


i*P 


P 

ac 


1 


where  P  is  the  orthogonal  projection  of  L2  (ft)  onto  K  . 
ac  sc 


§6.  The  Rayleigh-Bloch  Wave  Expansions  for  A. 

This  section  presents  a  construction,  based  on  the  results  of  §5, 
of  the  R-B  diffracted  plane  wave  eigenfunctions  ^+(X,p,q)  and  a  deriva¬ 
tion  of  the  corresponding  R-B  wave  expansions  for  A.  For  brevity  the 
derivation  is  restricted  to  the  cases  for  which  A  has  no  surface  waves; 
that  is,  a0(Ap)  *  <t>  for  all  p.  The  modifications  that  are  needed  when 
there  are  surface  waves  are  indicated  at  the  end  of  the  section. 

In  this  section  ^0<X,p,q)  denotes  the  R-B  wave  eigenfunction  for 
A0;  that  is,  one  of  the  functions  (1.33),  (1.34).  The  defining  proper¬ 
ties  of  t^+(X,p,q)  can  then  be  written 

(6.1)  6  D(A*0C)  ,  (p,q)  e  r*  , 

(6.2)  (A  +  U)2(p,q))  ^±(x,p,q)  ■  0  in  G  , 

(6.3)  <MX,p,q)  *  ^p(X,p,q)  +  ^(X,p,q)  in  R^  , 

where  (resp.,  ij/)  is  an  outgoing  (resp.,  incoming)  R-B  wave  for  G. 

The  construction  of  will  be  based  on  the  discussion  at  the 
end  of  §3.  Thus  If  (p,q)  €  R*  and  p  -  pQ  +  m  where  p0  €  (-1/2, 1/2]  and 
m  €  Z  then  ip+(X,p,q)  are  defined  by 

(6.4)  (X,p,q)  *  0P°<f>+(X,po+m,q)  , 
or,  more  explicitly, 

(6.5)  tp±(x,y,p,q)  *  exp  {27rUp0>  <t>+(x-27T£,y,p0+ra,q)  ,  (x,y)  €  . 

Theorem  5.1  then  implies 

59 

r  - - - - - 

J  E  d.  a  Cr  4  f  ^  Y  f  ,  L*y  ire 


►J*  •»*. 


Theorem  6.1.  Let  G  be  a  grating  domain  of  the  class  defined  in 
§1  and  let  A  »  A(G)  have  no  surface  waves.  Then  there  exist  unique  R-B 


diffracted  plane  waves  ^+(X,p,q)  for  each  (p,q)  e  R*  -  E,  where  E  is  the 
exceptional  set  (2.30).  Moreover,  the  mapping  (p,q)  +  ^+(*»P»q) 

€  ,^°C(A,G)  is  continuous  for  (p,q)  e  R*  -  E. 

The  principal  step  in  the  proof  of  Theorem  6.1  is  to  show  that 
defined  piece-wise  by  (6.5),  satisfies  (6.1).  This  may  be  done  by 
a  simple  distribution-theoretic  calculation  based  on  the  p-periodic 
boundary  condition  for  <J>+.  Details  are  given  in  [38].  The  uniqueness 
statement  follows  from  Theorem  2.1  since  o0(A^)  *  $  for  all  p  is 
assumed. 

The  R-B  wave  expansions  for  A  will  now  be  derived  from  the 
eigenfunction  expansions  for  A^  of  §5.  The  first  step  is  to  establish 
Parseval's  relation  for  A.  The  special  case  of  functions  f  E  L2°m(G)  is 
treated  first. 

Theorem  6.2.  For  all  f  e  L2°m(G)  define 


(6.6) 

f+(p*q)  * 

K(X.P.q)  f(X)  dX  ,  (p,q)  €  R?  -  E  . 

Then 

(6.7) 

f±  €  C(Rj  -  E)  n  L2 (Rj)  ,  and 

(6.8) 

' f ' L2 (G)  “  ,fz  L2(r£)  ‘ 

Proof . 

The  finiteness  of  f+(p,q)  for  (p,q)  €  R*  -  E  and  the 

A 

property  f±  E 

C(Rj  -  E) 

follow  from  the  last  statement  of  Theorem  6.1 

To  establish  the  rest  of  the  theorem  note  the  following  identity  for 


functions  f  E  L2°m(G)  and  points  (p,q)  E  R*  -  E. 


m 
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f+(p,q)  -  Mx.P.q)  f(x>  dx 
JG  " 

(6.9) 

=  l  I  m  K(x.p.q)  f(x)  dx 
stez  J£)w  “ 


“  l  <|>+  (x+2irJl,y,p,q)  f(x+2iri,,y)  dxdy 
SbeZ  JQ  ~ 

■  £  4>+(x,y,p,q)  e  2iriS,P  f(x+2u£.,y)  dxdy 

zez  Ja  - 

*  [  4>+(x,y,p,q)  [  l  e_27TiX,p  f (x+2ia,y)]  dxdy 


4>+(x,y,p,q)  F(x,y,p)  dxdy 

Jn  ~ 


where 


(6.10)  F(x,y,p)  -  l  e"27Tl^P  f(x+2*iT£,y)  ,  (x,y)  €  . 

j&sz 

Note  that  all  the  sums  in  (6.9)  are  finite  when  f  €  L,2°m(G),  Moreover, 

(6.10)  is  a  Fourier  series  in  p  with  a  fixed  finite  number  of  non-zero 
terms  for  all  (x,y)  €  Q. 

Equation  (6.9)  establishes  a  relation  between  the  eigenfunction 
expansions  for  A  and  Ap.  Thus  replacing  p  in  (6.9)  by  p  +  m  with 
p  €  (-1/2,1/21  and  m  €  Z  one  has 

(6.11)  f±(p+m,q)  -  F±(p+m,q,p) 

in  the  notation  of  §5.  In  particular,  Parseval's  relation  for  Ap, 
applied  to  F(* ,p) ,  and  (6.11)  gives 
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(6.12) 


|f(X,p) I2  dx 

a 


>00 


I 

m€Z  j0 


A 

f+(p-hn>q) 


2  dq 


Noting  the  continuity  of  p  F(*,p)  G  L2  (ft)  and  integrating  (6.12)  over 
p  €  (-1/2,1/21  gives 


f  1/2  t 

|F(X,p) I2  dx  dp  =  l 

■'-1/2  JQ  mez 


(6.13) 


1/2 


-1/2 


I  | f±(p+m,q) | 2  dqdp 
J  0 

L  |J±(p>q)|2  dpdq  =  "?±^2(r2)  • 


In  particular,  f+  €  L2(Rq)  which  completes  the  proof  of  (6.7).  To 
verify  (6.8)  note  that  Parseval's  formula  for  Fourier  series  implies 
that 


(6.14) 


l  |  f  (x+2irJl,y)  |  2  ,  X  €  Q  , 
2£Z 


where  the  sum  has  a  fixed  finite  number  of  terms  for  all  X  €  ft.  Inte¬ 
grating  (6.14)  over  X  €  (2  and  applying  Fub ini's  theorem  gives 


(6.15) 


[  [  |F(X,p)|2  dX  dp  -  l  f  |f  (x+2ir£,y)  j 2 

J -1/2  'ft  tez  'ft 

l  f  (SL)  If  (X)  |  2  dX 

ie z  Jn{'  ’ 


dX 


f  |£(X)|!  «  -  lfi^(G) 


Combining  (6.13)  and  (6.15)  gives  (6.8). 
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The  extension  of  Parseval's  relation  to  all  f  €  L2(G)  follows 
from  Theorem  6.2  by  a  standard  technique  using  the  denseness  of  L2°m(G) 
in  L2(G).  Thus,  writing 

(6.16)  (^  ■  G  n  (x  |  x2  +  y2  <  M2}  , 
one  has 

Corollary  6.3.  The  limits 

(6.17)  f  +  (p,q)  -  L2 (Rg)-lim  [  i|>+(X,p,q)  f(X)  dX 

M-«°  " 

exist  and  Parseval’s  relation  (6.8)  holds  for  all  f  €  1.2(G), 

A  representation  of  the  spectral  family  (II (p)  |  p  >  0}  of  the 
grating  propagator  A  will  now  be  derived  from  Corollary  6.3.  The  key 
fact  is  described  by 

Theorem  6.4.  The  resolvent  R(Afz)  *  (A  -  z)"1  of  the  grating 
propagator  A  satisfies  the  relation 

,  r  |f+(p.q)|2 

(6.18)  1  R(A,z)f | j  -  - = -  dpdq 

L*(G)  |u)Z(p,q)-z|2 

for  all  f  €  L2(G)  and  all  z  €  C  -  [0,®). 

To  prove  Theorem  6.4  it  is  enough  to  verify  (6.18)  for  all 
f  €  L2°m(G) .  The  idea  for  doing  this  is  to  define 

(6.19)  u(X)  -  R(A,z)  f(X) 

and  to  apply  Parsevalfs  relation  to  vM  «  <}>Mu  where  e  C^(R2).  For  a 
suitable  choice  of  one  has 
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(6.20)  vM  -  R(A,z) (f  +  gjj) 
where 

(6.21)  gjj  -  -2Vu  •  7^  ■  u  4^, 
and 

(6.22)  ^M±^P,q^  *  +  §h±(P»0)/w2(p»<I>  "  z 

whence 

(6.23)  l<^  R(A,z)fl  -  l(f±  +  g^/w2  ~  zl  • 

Passage  to  the  limit  M  ■+■  00  then  gives  (6.18).  For  the  case  of  the 
Dirichlet  boundary  condition  one  may  take  4^(X)  *  |x |  -  M)  where 
\l>  G  C°°(R)  satisfies  i//(x)  S  1  for  T  <  0  and  \p(x)  =  0  for  T  >  1.  For  the 
case  of  the  Neumann  boundary  condition  <f^  must  be  chosen  more  carefully, 
using  the  condition  G  €  S,  to  ensure  that  v^  satisfy  the  boundary 
condition.  The  details  of  the  construction  are  given  in  [38]. 

The  R-B  wave  expansions  for  A  follow  easily  from  Corollary  6.3 
and  Theorem  6.4.  They  are  formulated  as 

Theorem  6.5.  For  all  f  €  L2(G)  the  spectral  family 
(II (u)  |  U  >  0}  of  A  satisfies 


(6.24) 


n<p)  f(x) 


^+(X,p,q)  f±(p,q)  dpdq 


where 


Dy  -  R2  n  {(p,q)  |  p2  +  q2  <  y) 


(6.25) 


In  particular,  every  f  €  L2(G)  has  the  R-B  wave  expansion 


(6.26)  f (X)  =  L2 (G)-lim  ip+(X,p,q)  f+(p,q)  dpdq  . 

M-+*  J°M 

The  relation  (6.24)  is  a  direct  consequence  of  the  relation 

(6.27)  III(I)fl£  (G)  -  |  2  Xl^2(p,q))  |f±(p,q)|2  dpdq 

where  I  is  a  subinterval  of  [0,<»)  with  characteristic  function  Xj- 

(6.27)  follows  easily  from  (6.18)  and  Stone's  formula.  Note  that  (6.27) 
implies  the  absolute  continuity  of  the  grating  propagators. 

To  formulate  the  orthogonality  and  completeness  relations  for 
the  R-B  wave  expansions  define  linear  operators 

(6.28)  i±  :  L2(G)  -»•  L2(Rg) 
by 

(6.29)  $±  f  =  f±  • 

Then  and  are  spectral  mappings  for  A  in  the  sense  of 

Theorem  6.6.  For  every  bounded,  .Lebesgue-measurable  function 
¥(X)  defined  on  0  <  X  <  00 

(6.30)  $+  <F(A)  -  4((u)2(*))  $+ 

where  ¥(A)  is  defined  by  the  spectral  theorem. 

Moreover,  one  has 

Theorem  6.7.  The  R-B  wave  expansions  are  orthogonal  and  complete 
in  the  sense  that  and  are  unitary  operators: 
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(6.31)  *+  -  1  and  <I>+  $*  =  1 


Relations  (6.30)  and  the  completeness  relation  $+  $+  *  1  follow 
easily  from  the  spectral  theorem.  The  orthogonality  relation  $+  ®  1 

can  be  deduced  from  the  corresponding  property  of  $+p>  Theorem  3.8. 
Indeed,  it  is  sufficient  to  prove  that 


(6.32) 


<*+  *+f  -  f»f)  «  0 


for  all  f  in  a  dense  subset  of  L2(Rq).  This  may  be  verified  by  direct 
calculation  using  f  €  C0(R<)  -  E)  and  the  orthogonality  relation  for 
The  details  are  given  in  [38], 

Operators  A  that  admit  R-B  Surface  Waves.  It  was  shown  in  §2 

that  for  each  p  6  (-1/2, 1/2]  A  may  have  R-B  surface  waves  ^(X,p)  and 

eigenvalues  A.(p)  with  x-momentum  p.  The  functions  4>.(X,p)  =  ^.(X,p)[Q 
J  J  J  M 

are  precisely  the  eigenfunctions  of  A^.  The  principal  difficulty  in 
constructing  an  eigenfunction  expansion  for  A  in  this  case  is  in  con¬ 
structing  families  of  R-B  surface  waves  ^ (X,p)  and  eigenvalues  A^ (p) 
whose  dependence  on  p  is  sufficiently  regular.  The  "axiom  of  choice11 
definition  (independent  choice  for  each  p)  in  inadequate  to  give  even 
measurability  in  p.  This  was  pointed  out  in  the  author’s  paper  on  the 
analogous,  but  simpler,  case  of  Bloch  waves  in  crystals  [36]. 

If  9 G  is  a  union  of  smooth  curves  (class  C3)  then  the  Green’s 
functions  (4.44),  (4.45)  can  be  used  to  construct  an  integral  equation 
for  the  eigenfunctions  <}>j(X,p).  In  this  case  the  method  of  [36]  can  be 
used  to  construct  ’’almost  holotnorphic”  families  (<t>j(X,p)}. 

In  the  general  case  there  is  a  one-to-one  correspondence  between 

eigenfunctions  <p.  (X,p)  of  A  and  eigenfunctions  9,(X,p)  of  A  with 

J  P  J  P»s*r 
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eigenvalues  n  (?)  €  [p2,°°)  given  by  0.(*,p)  =  P  .  <  * ,  p)  .  The 

r  J  P  ><■■>>  r  J 

eigenvalues  of  ^  ^  are  isolated,  with  finite  multiplicity,  and  may 
be  studied  by  the  methods  of  analytic  perturbation  theory  (Kato  [17, 

Ch .  7]).  These  problems  will  not  be  pursued  here. 

If  a  sufficiently  regular  family  of  R-B  surface  waves  for  A  has 
been  constructed  the  eigenfunction  expansions  for  A  may  be  derived  by  the 
method  introduced  above.  Thus,  defining  ^(X,p)  =  0  when  j  >  N(p)  , 
equation  (6.12)  must  be  replaced  by 


(6.33) 


If(X.p) l 2  dx  =  l  |f  (P) [ 2  +  l 

j*l  J  m€Z  * 


|f+(p+m,q) | 2  dq 


where 


(6.34) 


VP) 


iMx.p)  f(x)  dx  • 


Integration  over  p  6  (-1/2, 1/2]  gives  the  Parseval  relation 


(6.35) 


,f,U(G) 


£  'VV-I/!,!/!)  +  ' 


The  corresponding  representation  of  the  spectral  family  is 


n(u)  f(x) 


r 1/2 

l  <l>.(X,p)  f.(p)  dp 
-1/2  (p)<y  J  J 


+ 


^+(X,p,q)  f±(p,q)  dpdq  . 


(6.36) 


§7.  Concluding  Remarks. 


It  is  clear  that  the  methods  developed  in  this  report  are 
applicable  to  diffraction  gratings  in  R3  (and  in  Rn,  n  >  3) .  They  are 
also  applicable  to  gratings  with  holes,  as  in  Alber  [3]  where  the 
grating  Rn  -  G  is  a  periodic  structure  contained  in  a  slab  R^  ,  and 
G  is  connected.  The  changes  needed  to  treat  these  cases  are  primarily 
notational.  The  same  methods  can  also  be  applied  to  the  physically 
important  cases  of  dielectric  gratings  in  electromagnetic  theory  and 
elastic  gratings  in  acoustics. 

The  most  important  unsolved  problems  for  diffraction  gratings 
concern  the  R-B  surface  waves.  Gratings  that  admit  such  surface  waves 
would  presumably  be  good  waveguides  for  signals  generated  near  the 
grating  surface.  Geometric  criteria  for  the  existence  and  non-existence 
of  such  waves  would  be  of  great  interest  for  applications. 
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